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Yield  Tables  for  Managed  Even- Aged  Stands 
of  Spruce-Fir  in  the  Central  Rocky  Mountains 


Robert  R.  Alexander,  Wayne  D.  Shepperd, 
and  Carleton  B.  Edminster 


Timber  management  in  the  Engelmann 
spruce  {Picea  engelmannii)-suba\pine  fir 
(Abies  lasiocarpa)  type  in  the  central  and 
southern  Rocky  Mountains  is  in  a  period  of 
transition  to  more  intensive  and  varied 
management  practices.  Large  areas  of  old- 
growth  forests  are  rapidly  being  converted  to 
stands  that  must  be  managed  from  the 
regeneration  period  to  harvesting.  Yield 
tables  that  predict  probable  yields  of  wood 
that  will  result  from  specified  combinations 
of  site  quality,  frequency  and  intensity  of 
thinning,  and  utilization  standards  provide 
goals  toward  which  conversion  can  be 
directed.  They  also  provide  part  of  the  infor- 
mation needed  to  determine  the  influence  of 
timber  management  practices  on  other  forest 
resources  (Myers  1971).  Land  managers  can 
examine  alternatives  and  make  decisions 
through  computerized  evaluation  (Program 
TEVAP2)  of  forestry  activities,  a  procedure 
that  uses  yield-table  computation  in  a  set  of 
mathematical  operations  (Myers  1974). 

Procedures  for  the  computation  of  yield 
tables  for  managed  even-aged  stands  of 
spruce-fir  presented  in  this  Paper  were  ad- 
apted from  the  field  and  computer  procedures 
for  managed  stand  yield  tables  (Program 
PONYLD)  developed  by  Myers  (1971).  In- 
cluded are  (1)  identification  of  the  field 
measurements  that  provide  the  basic  data 
needed  to  produce  the  yield  tables,  (2)  the 
relationships  that  replace  the  species-specific 
statements  in  Myers'  Program  PONYLD,  (3) 
the  computer  program  (SPRYLD)  written  in 
FORTRAN  IV  that  computes  and  prints 
yield  tables  for  spruce-fir,  and  (4)  an  example 
of  what  program  SPRYLD  can  produce.  The 
computer  program  SPRYLD  has  the  capacity 
of  producing  a  series  of  yield  tables  which 
show  how  projected  outcomes  will  vary  in 
response  to  changes  in  cultural  treatments 
and/or  variations  in  original  stand  and  site 


conditions.  Large  numbers  of  tables,  each 
based  on  a  specific  set  of  alternatives,  can  be 
computed  and  printed  for  a  few  cents  each. 
They  provide  the  manager  with  the  oppor- 
tunity to  examine  the  probable  results  of  his 
operations,  make  necessary  changes  in 
management  goals,  and  study  the  effect  of 
these  changes  before  money  is  spent  on  them 
(Myers  1971). 

General  Description  of  Methods 

The  nine  items  of  information  (working 
tools)  described  later  are  species-specific 
statements  in  Myers'  (1971)  computer 
program.  All  but  the  first  item,  the  desired 
residual  stand  after  each  cutting  which  is 
based  on  available  information  from  thinned 
stands,  were  replaced  to  compute  yield  tables 
for  even-aged  managed  stands  of  spruce-fir. 
Data  for  these  eight  items  of  information 
were  based  on  temporary  growth  prediction 
plots  measured  in  detail.  Tree  volume 
equations  (Myers  and  Edminster  1972)  and 
site  index  curves  (Alexander  1967)  applicable 
to  spruce-fir  stands  in  the  central  Rocky 
Mountains  were  used  in  the  analyses  of  data 
described  in  the  following  sections. 

There  are  no  thinned  stands  of  spruce-fir 
in  the  central  Rocky  Mountains,  and  partial- 
ly cut  stands  were  unsuitable  for  sampling; 
therefore,  the  growth  prediction  plots  were 
placed  in  69  undisturbed  even-aged  stands 
throughout  the  spruce-fir  type  in  Colorado 
and  Wyoming.  These  plots,  chosen  to  ap- 
proximate what  a  managed  stand  might  look 
like,  conformed  to  the  usual  requirements  as 
to  uniformity  of  site  quality,  range  in  tree 
sizes,  and  stand  density  across  the  plot.  The 
trees  were  also  free  of  diseases  or  insect 
infestations  that  would  affect  growth.  The 
plots  covered  a  range  of  site  quality  (SI  40  to 


1 


110)  and  stand  density  (30  to  1,600  stems  per 
acre).  Age  classes  varied  from  20  to  160  years. 

Measurements  made  on  each  growth 
prediction  plot  included  the  following: 

1.  Plot  area. 

2.  Heights  and  ages  of  6  to  10  dominant  trees 
suitable  for  site  index  determination 
(Alexander  1967). 

3.  Diameter  at  breast  height  (d.b.h.)  of  each 
tree  to  the  nearest  0.1  inch. 

4.  Total  height  of  each  tree  to  the  nearest  1.0 
foot,  or  a  sufficient  sample  in  each 
diameter  class  to  construct  a 
height/ diameter  curve  where  a  large 
number  of  trees  were  measured.  Heights  of 
all  dominants  and  codominants  were 
measured  except  where  tops  were  dead, 
defective,  or  deformed. 

5.  Crown  class  of  each  tree. 

6.  Total  ages  of  a  sample  of  dominant  and 
codominant  trees  to  validate  even-age 
status  of  the  main  stand. 

7.  Radial  wood  growth  during  the  past  10 
years  from  increment  borings  at  breast 
height  along  an  average  radius. 

8.  D.b.h.  outside  bark  of  all  trees  that 
appeared  to  have  died  during  the  last  10 
years. 

9.  On  12  plots,  a  number  of  cut  or  leave  codes 
for  each  tree  based  on  trial  marking  to 
simulate  several  intensities  of  thinning. 
This  information  was  used  to  provide 
some  of  the  input  needed  to  determine 
diameter  and  height  increases  due  to  thin- 
ning, stand  volumes,  and  volume  conver- 
sion factors.  Additional  simulated  data 
were  generated  using  Program  SPRCHK, 
a  modification  of  Program  PONCHK 
(Myers  1971).  The  12  plots  were  also  used 
to  verify  the  SPRCHK  output. 

Plot  and  tree  data  computed  initially 
from  the  field  measurements  included: 

1.  Site  index  (Alexander  1967). 

2.  A  height-diameter  curve  for  each  plot  to 
provide  a  height  for  each  tree  for  which 
actual  height  was  not  measured. 

Measured  and  computed  items  that 
describe  the  present  stand  were  used  to 
compute  the  following  values  for  each  plot: 

1.  Number  of  trees  per  acre. 

2.  Number  of  dominants  and  codominants 
per  acre. 

3.  Basal  area  per  acre. 

4.  Average  d.b.h.,  computed  as  the  tree  of 
average  basal  area. 

5.  Average  height  of  dominant  and  codomi- 
nant trees. 


6.  Average  height  of  all  trees. 

7.  Average  main  stand  age. 

8.  Total  cubic  feet  from  ground  to  tip  for  all 
trees,  per  acre. 

9.  Merchantable  cubic  feet  to  a  4-inch  top  in 
trees  5  inches  d.b.h.  and  larger,  per  acre. 

10.  Board  feet  (Scribner  Rule)  to  a  6-inch  top 
in  trees  8  inches  d.b.h.  and  larger,  per 
acre. 

All  except  item  7  were  generated  as  part 
of  the  output  of  Program  SPRCHK. 

Diameters  of  live  trees,  diameters  of  the 
tallied  dead  trees,  and  present  stand  age  on 
each  plot  provide  the  following  items  that 
described  the  stands  10  years  ago  at  the 
beginning  of  the  prediction  period. 

1.  Number  of  trees,  per  acre. 

2.  Basal  area,  per  acre. 

3.  Past  d.b.h.  of  each  tree  from  present  d.b.h., 
radial  wood  growth,  and  periodic  bark 
growth  (Myers  and  Alexander  1 972). 

4.  Average  stand  diameter  (tree  of  average 
basal  area). 

5.  Average  main  stand  age  (present  main 
stand  age  minus  10  years). 

Development  of  Items  to  Replace 
Species-Specific  Statements  and 
Other  Modifications  Needed  to 
Adapt  Myers'  (1971)  Program 
PONYLD  to  SPRYLD 

After  plot  measurements  were  obtained 
and  summarized,  the  items  described  below 
were  computed  as  one  or  more  relationships 
to  convert  the  species-specific  statements  in 
Program  PONYLD  to  spruce-fir.  Most  of  the 
nine  items  that  appear  as  FORTRAN 
statements  in  Program  SPRYLD  and  its 
associated  subroutines  SPRCUT  and  SPR- 
VOL  were  obtained  by  regression  analysis  of 
plot  values  described  above. 


1.   Stocking  After  Cutting 

Stand  density  to  be  left  after  each  cutting 
is  expressed  as  the  relationship  between 
basal  area  and  average  stand  diameter.  Data 
on  spruce-fir  from  thinning  studies  or  tem- 
porary plots  were  not  available  to  construct  a 
graph  of  desired  basal  area  over  stand 
diameter  for  local  average  diameter.  Data  for 
this  item  were  taken  directly  from  Myers' 
(1971)  Program  PONYLD.  When  sufficient 


information  for  spruce-fir  becomes  available, 
this  item  will  be  examined  and  necessary 
adjustments  made. 

The  following,  taken  from  Myers  (1971), 
is  repeated  here  for  continuity. 

"In  table  1,  basal  area  increases  with 
diameter  until  10.0  inches  diameter  is  reach- 
ed, and  remains  constant  thereafter.  The 
designation  'growing  stock  level  8  0'  indicates 
that  basal  area  is  80.0  ft2  when  diameter  is 
10.0  inches  or  larger,  regardless  of  what  basal 
area  may  be  at  lower  average  diameters. 
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Figure  1.  —  Basal  area  after  thinning  in 
relation  to  average  stand  diameter  for 
standard  levels  of  growing  stock  (from 
Myers  1971). 


"Desired  stand  density  will  vary  with  the 
objectives  of  management,  and  a  family  of 
basal  area-diameter  relationships  is  needed 
(fig.  1).  The  original  single  curve  or  function 
of  basal  area  on  diameter  is  treated  as  a  guide 
curve  from  which  other  curves  can  be  produc- 
ed. Basal  areas  for  any  growing  stock  level 
can  be  computed  by  multiplying  the  guiding 
values  for  level  80  in  table  1  by  the  ratio 
level/80." 

The  level  designations  that  are  the 
variables  THIN,  REST,  and  DSTY  in 
SPRYLD  are  the  same  as  appear  in 
PONYLD. 

The  curves  of  figure  1  define  growing 
stock  goals  for  many  possible  management 
objectives.  Any  desired  form  of  the  guide 
curve  may   be   used  if  the  appropriate 


Table  1. --Basal  areas  after  intermediate  cutting 
in  relation  to  average  stand  diameter  growing 
stock  level  80  (from  Myers  1971) 
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statements  of  subroutine  SPRCUT  are 
modified  properly. 

Relationships  shown  in  table  1  appear  as 
functions  for  level  80  in  program  SPRYLD. 
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Basal  areas  computed  from  these  functions 
are  multiplied  by  terms  that  include  the 
desired  growing  stock  level  (THIN)  to  obtain 
values  for  other  growing  stock  levels. 
Variables  for  which  FORTRAN  statements 
were  taken  from  PONYLD  and  their  use, 
were: 

"a.  DBHP  -  to  find  a  d.b.h.  less  than  10.0 
inches  when  basal  area  is  known.  Three 
equations  for  DBHP  are  used  to  simplify 
representation  of  the  nonlinear 
relationship  between  d.b.h.  and  basal 
area. 

b.  BREAK  and  BUST  —  to  compute  values 
of  basal  area  that  are  the  upper  limits  of 
applicability  of  the  first  two  equations  for 
DBHP. 

c.  SQFT  —  to  find  basal  area  when  d.b.h.  is 
known.  Two  equations  represent  the  non- 
linear relationship  for  d.b.h.  less  than 
10.0  inches." 


with  an  average  diameter  of  4.4  inches  for  site 
index  classes  40  and  50.  As  more  young 
stands  with  reasonable  spacing  reach  30 
years  of  age,  additional  data  will  become 
available  to  evaluate  this  first  approxima- 
tion, and  make  necessary  adjustments. 


3.  Diameter  Increase  from  Growth 

Regression  analysis  of  data  from  growth 
prediction  plots  provided  an  equation  for 
predicting  future  average  stand  diameter  of 
spruce-fir.  Present  average  stand  diameter  is 
estimated  from  past  average  stand  diameter, 
site  index,  and  past  basal  area  per  acre.  The 
following  equation  for  a  10-year  prediction 
period  appears  in  SPRYLD  as  the  FORTRAN 
statement  for  DBHO: 


2.  Description  of  Unthinned,  Young 
Stands 

Values  in  the  first  line  of  each  yield  table 
describe  stand  conditions  just  prior  to  initial 
thinning.  They  are  entered  directly  from  data 
cards  or  are  computed  from  the  data.  Users  of 
SPRYLD  must,  therefore,  be  able  to  describe 
the  stands  that  do  or  should  exist  at  time  of 
initial  thinning.  In  a  yield  table  for  managed 
stands,  the  stand  density  and  related  average 
diameter  given  in  the  first  line  result  when 
stand  regeneration  and  subsequent  growth 
and  mortality  progress  as  planned. 

Only  a  few  unthinned  young  spruce-fir 
stands  were  found  that  represented  possible 
regeneration  goals  for  various  management 
objectives.  Furthermore,  no  usable  data  were 
available  from  yield  studies  made  elsewhere. 
It  was  not  possible,  therefore,  to  determine 
average  diameter  for  each  site  class  for 
various  combinations  of  age  and  number  of 
trees  per  acre.  Instead,  an  average  diameter 
in  relation  to  number  of  trees  per  acre,  was 
determined  for  an  average  site  index  at  age  30 
years,  the  youngest  age  appearing  in  the  yield 
tables.  It  was  then  necessary  to  make  a 
judgment  decision  as  to  the  average 
diameters  and  number  of  trees  that  appeared 
reasonable  for  managed  stands  at  the  time  of 
first  thinning  for  each  site  class.  These  values 
range  from  800  stems  per  acre  with  an 
average  diameter  of  4.7  inches  for  site  index 
classes  100  and  greater  to  950  stems  per  acre 


DBHO  =  1.62917  +  1.03371  (DBHT) 
+  0.01304  (SITE)  -  0.90669  (log10  BAST) 
Syx  =  0.2469   R2  =  0.9948 

where      DBHO  =  present  average  stand  dia- 
meter. 

DBHT  =  past  average  stand  dia- 
meter. 
SITE  =  site  index. 
BAST  =  past  basal  area  per  acre. 


4.  Diameter  Increase  from  Thinning 

The  change  in  average  stand  diameter 
resulting  from  intermediate  cuttings  was  es- 
timated from  data  obtained  by  trial  marking 
of  growth  prediction  plots,  and  the 
supplemental  procedure  developed  by  Myers 
(1971)  and  generated  by  Program  SPRCHK. 
This  later  procedure  provided  simulated  data 
on  combinations  of  initial  stand  diameters, 
stocking  level,  and  stand  density  not 
available  from  the  trial  markings. 

In  subroutine  SPRCUT,  diameter  after 
thinning  is  estimated  from  diameter  before 
thinning  and  the  percentage  of  trees  to  be 
retained.  Regression  analysis  of  data  from 
simulated  thinnings  provides  two  functions 
—  DBHE  and  PDBHE  —  that  represent  the 
same   variable,   diameter  after  thinning. 
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DBHE  is  computed  directly  in  subroutine 
SPRCUT  from  the  following  equation  if  at 
least  50  percent  of  the  trees  are  retained: 

DBHE  =  0.02666  +  1.30655  (DBHO)  - 
0.00306  (DBHO  x  PRET) 


yx 


0.1413,  R2  =  0.9985 


where      DBHE  =  average    diameter  after 
thinning. 
DBHO  =  average  stand  diameter 

before  thinning. 
PRET  =  percentage   of  trees  re- 
tained after  thinning. 

With  fewer  trees  retained,  the  relationship  is 
nonlinear,  so  PDBHE  is  computed  in  sub- 
routine SPRCUT  from  the  following  equation 
and  its  antilogarithm  becomes  DBHE: 

PDBHE  =  0.33206  +  0.98346  (log10  DBHO) 
-  0.14170  (log10  PRET) 

S„Y  =  0.0187,   R2  =  0.9898 
yx 


5.  Heights  of  Dominants  and 
Codominants 

Average  heights  of  dominant  and 
codominant  trees  were  computed  from  data 
from  the  growth  prediction  plots  and  from 
Alexander's  (1967)  site  index  curves  adjusted 
from  dominant  height  to  dominant  and 
codominant  height  (table  2).  Regression 
analysis  provided  two  equations  that  appear 
as  statements  for  HTSO  in  SPRYLD,  for 
estimating  height  for  various  combinations 
of  age  and  site  index: 

The  equations  are  shown  below: 

HTSO  (AGEO  <  100)  =  -  13.71751  + 
0.15087  (SITE)  +  0.00126  (AGEO2)  + 
0.01371  (AGEO  x  SITE) 
-  0.00006  (AGEO2  x  SITE) 

SVY  =  2.0074,    R2  =  0.9920 
yx 

log10  HTSO  (AGEO  >  100)  =  0.91859 
-  10 AGEO-1  +  062318  (logl°  SITE) 


Table  2. — Average  height  (feet)  of  dominant  and 
codominant  trees  at  various  ages  and  site 
indexes  for  Engelmann  spruce  (as  computed  by 
statements  for  HTSO) 


Ma  i  n 

Site 

index 

(Based  on 

domi  nant 

trees) 

stand 

age 

(years) 

50 

60 

70 

80 

90 

20 

3 

7 

1  1 

15 

19 

23 

27 

30 

8 

13 

1  o 
1  0 

23 

2o 

33 

38 

4o 

12 

18 

25 

31 

37 

43 

49 

50 

17 

24 

31 

37 

1. 1. 

44 

51 

58 

60 

21 

29 

36 

hh 

51 

59 

67 

70 

25 

33 

Hi 

50 

58 

66 

74 

80 

29 

37 

It  L 

46 

55 

63 

72 

81 

90 

33 

41 

50 

60 

68 

77 

87 

100 

36 

45 

54 

64 

73 

82 

92 

l  i  n 

3? 

10 

58 

67 

77 

86 

120 

41 

51 

61 

70 

80 

89 

99 

130 

43 

54 

63 

73 

83 

92 

102 

HO 

45 

56 

66 

76 

86 

95 

105 

150 

47 

58 

68 

78 

88 

98 

107 

160 

49 

60 

70 

80 

90 

100 

109 

170 

51 

61 

72 

82 

92 

101 

111 

180 

53 

63 

73 

83 

93 

103 

113 

Syx  =  °-0101' 


R2  =  0.9938 


where       HTSO  =  average  height  of  domi- 
nant   and  codominant 
trees  before  thinning. 
SITE  =  site  index. 
AGEO  =  main  stand  age. 


6.  Height  Increase  from  Thinning 

Changes  in  the  average  height  of  domi- 
nant and  codominant  trees  due  to  in- 
termediate cuttings  were  estimated  from  data 
provided  by  trial  markings  of  growth  predic- 
tion plots,  and  simulated  thinnings 
generated  by  Program  SPRCHK  for  various 
additional  combinations  of  stocking  level, 
initial  stand  diameter,  and  initial  stand  den- 
sity. In  SPRYLD,  the  variable  ADDHT  is  the 
computed  amount  of  change  estimated  from 
the  percentage  of  trees  retained  (PRET)  by 
the  equation  below: 


+ 


40.08154  logX(fEOE 


ADDHT  =  -1.81392  +  17.77922.rZi_ 

'V  PRET 

Sv  =  0.0864,   R2  =  0.9958 
yx 
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At  each  cutting,  the  current  value  of 
ADDHT  is  added  to  height  before  thinning, 
HTSO,  to  obtain  height  after  thinning, 
HTST.  It  is  also  added  to  a  cumulative  sum  of 
changes,  HTCUM,  so  that  computed  heights 
before  thinning  will  show  the  effects  of  past 
treatments  as  well  as  of  age.  Values  for 
ADDHT  are  small  because  changes  in  height 
are  for  dominant  and  codominant  trees  only. 


7.  Noncatastrophic  Mortality 

Mortality  in  unthinned  stands  is  usually 
more  important  than  in  thinned  stands.  Since 
the  only  data  available  were  from  growth 
prediction  plots,  and  they  were  located  in 
unthinned  stands  with  spacings  that 
appeared  desirable  for  managed  stands,  es- 
timates of  mortality  in  managed  stands  are 
only  approximate.  Furthermore,  we  did  not 
compute  a  prediction  equation  for  mortality 
for  spruce-fir  stands  with  an  average 
diameter  of  10.0  inches  or  larger  because  of 
the  wide  variability  in  mortality  in  those 
stands. 

The  prediction  equation  in  SPRYLD  for 
stands  with  an  average  diameter  of  less  than 
10.0  inches,  shown  below  for  the  percentage 
of  mortality  expressed  as  a  decimal  (DIED), 
contains  average  stand  diameter  and  basal 
area  at  the  beginning  of  the  period  as  in- 
dependent variables: 

DIED  =  -  0.0003967  + 

0.0000382  (DBHT  x  BAST) 


Sv  =  0.0055,   R2  =  0.9677 


simulated  thinnings,  and  trial  markings  were 
converted  to  total  cubic  foot,  merchantable 
cubic  foot,  and  board  foot  volumes  per  acre 
with  appropriate  tree  volume  equations 
(Myers  and  Edminster  1972).  In  subroutine 
SPRVOL,  only  total  cubic  foot  volumes  per 
acre  are  computed  directly,  and  therefore  are 
the  only  unit  of  measure  for  which  a  stand 
volume  equation  is  needed. 

Stand  volume  equations  appear  as 
statements  to  compute  CUFT  in  subroutine 
SPRVOL.  Two  statements  were  used  because 
the  relationship  was  not  linear  over  the  range 
of  D2H  that  can  appear  in  the  yield  tables: 


CUFT  (D2H  <  22,500)  =  (0.344  B0  - 
0.00524BA  +  0.0023575D2H)  x  N 

Svv  =  0.4982,   R2  =  0.99  86 
yx 

CUFT  (D2H  >  22,500)  =  (7.60196  - 
0.01052BA  +  0.0020103D2H)  x  N 

Syx  =  1.9210,   R2  =  0.9940 

where        CUFT  =total  cubic  foot  volume  per 
acre. 

D=  average  stand  diameter. 
BA=  basal  area  per  acre. 

H=  average  height  of  domi- 
nant and  codominant 
trees. 

N=  number  of  trees  per  acre. 


The  standard  errors  of  estimate  quoted  above 
result  from  the  prediction  of  cubic  foot  volume 
of  a  tree  of  average  d.b.h.  and  height,  before 
multiplying  by  number  of  trees  per  acre. 


where  DIED  =  percentage  of  trees  that 
have  died  during  mea- 
surement period  (10  years). 
DBHT  =  average  stand  diameter  at 
the  beginning  of  the  mea- 
surement period. 
BAST  =  basal  area  per  acre  at  the 
beginning  of  the  measure- 
ment period. 

8.  Stand  Volume  Equation 

Plot  tallies  of  tree  diameters  and  heights 
from  growth  prediction  plots,  the  result  of 


9.  Volume  Conversion  Factors 

Subroutine  SPRVOL  computes  volumes 
in  merchantable  cubic  feet  and  board  feet 
from  total  cubic  volume  per  acre  and  ap- 
propriate conversion  factors.  Data  from 
growth  prediction  plots,  trial  markings,  and 
simulated  thinnings  that  produced  the  stand 
volume  equations  also  provided  conversion 
factors.  The  quantity  of  each  unit  per  total 
cubic  foot  was  determined  separately  for  each 
data  point.  A  stand  diameter  of  5.0  inches  and 
top  diameter  of  4.0  inches  for  merchantable 
cubic  feet,  and  a  stand  diameter  of  8.0  inches 
and  top  diameter  of  6.0  inches  for  board  feet, 
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were  selected  as  appropriate  minimums. 
Regression  analysis  provided  the  function  to 
compute  the  conversion  from  total  cubic 
volume  to  merchantable  cubic  volume  and 
board  foot  volumes  that  appear  as  statements 
FCTR  and  PROD  in  SPRVOL.  FCTR  can  be 
estimated  by  the  following  equation  from 
average  stand  diameter: 


FCTR  =  0.82375  + 


3.45569 


+  0.00013(D2)  - 


D 

28.86783 

Do 


Syx  =  0.0273,   R2  =  0.9572 


Estimates  of  PROD  for  spruce-fir  can  be 
improved  if  basal  area  is  included  with 
diameter  in  the  equations.  Two  equations 
shown  below  appear  in  SPRVOL  so  the 
relationships  can  be  expressed  in  simpler 
terms  over  a  wide  range  of  diameters: 


PROD  (D  <16.5  inches)  =  4.59159  - 

214-06370  +  0.39782  (log10  BA) 

S     =  0.1567,   R2  =  0.9561 
yx 

PROD  (D  >16.5  inches)  =  8.59422  - 

19-54507  +  0.44432  (log10  BA) 
vD 

Syx  =  0.1191,   R2  =  0.8575 

Description  of  Program  SPRYLD 

Program  SPRYLD  consists  of  a  main 
program  and  two  subroutine  subprograms 
written  in  standard  FORTRAN  IV.  The  main 
program  reads  the  data  cards,  performs  most 
computations,  and  writes  the  yield  tables. 
Subroutine  SPRCUT  determines  the  new 
average  stand  diameter  and  percentage  of 
trees  retained  after  cutting  to  the  specified 
growing  stock  level.  Subroutine  SPRVOL 
computes  volumes  in  total  cubic  feet  per  acre, 
and  factors  to  convert  these  volumes  to  other 
units.  Operations  performed  by  each  routine 
are  indicated  by  comment  statements  in  the 
source  listing  (Appendix  1). 

As  mentioned  earlier,  program  SPRYLD 
computes  and  prints  sets  of  yield  tables.  The 


variable  NTSTS,  read  from  data  card  type  1, 
controls  the  number  of  sets  of  tables.  The 
number  of  yield  tables  within  a  set  is  con- 
trolled by  variable  MIX,  read  from  data  card 
type  3.  The  first  yield  table  of  each  set  is 
computed  from  initial  conditions  and  con- 
trols on  operations  specified  on  data  cards 
type  3,  4,  and  5.  Subsequent  yield  tables 
within  a  set  are  computed  from  the  same 
conditions  and  controls,  with  the  exception  of 
the  growing  stock  level  for  intermediate  cuts 
after  initial  thinning.  This  growing  stock 
level  is  increased  by  the  value  of  DSTINC 
(read  from  data  card  type  3)  over  the  level  of 
the  previous  yield  table.  Operations  per- 
formed for  each  yield  table  are: 

1.  Computation  of  basal  area,  height,  and 
volumes  just  prior  to  initial  thinning. 

2.  Partial  cutting  to  the  growing  stock  level 
specified  for  initial  thinning  or  subsequent 
cutting.  Cutting  will  not  be  simulated  if  the 
stand  is  already  below  the  growing  stock 
level  specified. 

3.  Computation  of  post-cutting  density,  basal 
area,  height,  and  volumes. 

4.  Printing  of  table  headings  once  for  each 
yield  table  and  printing  of  values  ap- 
propriate for  the  stand  age. 

5.  Projection  of  diameter,  height,  and  stand 
density  for  one  or  more  periods  until  the 
next  intermediate  cut  is  scheduled.  Stand 
volumes  and  other  values  are  computed 
and  printed  at  ages  when  no  cutting  is 
scheduled. 

6.  Repetition  of  steps  2  to  5  until  stand  age  at 
time  of  initial  regeneration  cutting  is 
reached. 

7.  Redefinition  of  the  growing  stock  to  be  left 
after  cutting  and  the  interval  between 
regeneration  cuttings  if  shelterwood  cut- 
tings are  specified.  Table  computations 
will  be  terminated  if  clearcutting  is 
specified. 

8.  Repetition  of  steps  2  to  7  to  accomplish 
regeneration  cuts  until  the  age  of  the  final 
cut  is  reached. 

9.  Printing  of  totals  for  volumes  removed. 
Volumes  less  than  COMCU  or  COMBF 
(read  from  data  card  type  2)  will  not  be 
included  in  the  totals  so  that  total  commer- 
cial yields  may  be  compared,  if  desired. 
Actual  totals  may  be  obtained  by  entering 
values  of  1.0  for  COMCU  and  COMBF  on 
data  card  type  2., 

Subroutine  SPRCUT  computes  average 
stand  diameter  after  cutting  from  diameter 
before  cutting  and  the  percentage  of  trees 
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retained.  The  percentage  of  trees  retained  is 
needed  as  an  independent  variable,  but  is 
itself  an  unknown.  Successive  percentages  of 
trees  retained  are  tested  until  d.b.h.  after 
thinning,  number  of  trees  retained,  and 
residual  basal  area  agree  with  the  diameter 
and  basal  area  combination  specified  by  the 
growing  stock  level  of  the  cut.  These  com- 
binations are  shown  in  table  1  for  growing 
stock  level  80.0. 

Subroutine  SPRVOL  computes  total 
cubic  feet  per  acre  and  factors  to  convert  this 
volume  to  other  units.  Conversions  to 
merchantable  cubic  feet  and  to  board  feet  are 
shown  in  Appendix  1.  Utilization  standards 
for  these  units  are  given  in  subsection  9  of  the 
previous  section  and  in  the  comment 
statements  of  SPRVOL.  Conversions  to  other 
units  or  utilization  standards  may  supple- 
ment or  replace  those  already  in  SPRVOL. 

Program  SPRYLD  should  run  with  little 
or  no  modification  on  any  computer  that 
accepts  FORTRAN  IV,  has  a  minimum  of 
32K  words  of  storage,  and  has  two  input/out- 
put devices  (unit  5  for  program  and  data  deck 
input  and  unit  6  for  printed  output).  Changes 
to  adapt  the  program  to  other  utilization 
standards  and  for  additional  computations 
are  described  by  Myers  (1971)  in  his  section 
headed  Modifications  of  PONYLD,  and  are 
not  repeated  here. 


Description  of  the  Data  Deck  for 
Program  SPRYLD 

The  data  deck  for  program  SPRYLD 
consists  of  five  types  of  data  cards.  These 
cards  are  numbered  by  their  order  of 
appearance  in  the  data  deck  except  that 
cards  type  3,  4,  and  5  may  be  repeated  in  sets 
as  specified  by  variable  NTSTS  from  card 
type  1.  The  first  two  cards  of  the  data  deck 
(types  1  and  2)  enter  values  which  do  not 
change  during  a  computer  run.  Cards  type  3, 
4,  and  5  enter  values  used  in  the  computation 
of  a  set  of  yield  tables  which  may  change 
between  sets  of  tables.  The  contents  of  each 
data  card  are  described  in  the  following 
tabulation  of  the  order  and  contents  of  the 
data  deck. 


SPRYLD  and  the  printed  results  obtained.  It 
illustrates  some  of  the  questions  that  may  be 
asked  and  the  information  that  will  be 
provided.  The  example  also  serves  as  a  test 
problem  for  use  in  adapting  the  source 
program  to  locally  available  computing 
facilities. 

A  forest  manager  wishes  to  determine 
the  intensity  of  thinning  that  will  maximize 
volume  production  in  board  feet  in  stands  of 
site  index  80  (SITE).  Length  of  the  cutting 
cycle  (JCYCL)  has  not  been  standardized, 
but  will  be  30  years  for  this  test.  He  also 
wants  to  compare  yields  from  two-cut  and 
three-cut  shelterwood,  both  with  the  final 
removal  cut  scheduled  for  stand  age  150 
years  (REGN(2)  for  two-cut  shelterwood, 
REGN(3)  for  three-cut  shelterwood)  and  con- 
sidering the  current  crop  only.  Alternatives 
calling  for  more  than  one  precommercial 
thinning  are  unacceptable.  Minimum  com- 
mercial volumes  per  acre  are  400  cubic  feet 
(COMCU)  and  2000  board  feet  (COMBF).  The 
manager  expects  that  his  procedure  for 
regeneration  cuts  will  result  in  a  new  stand 
that  contains  850  trees  per  acre  (DENO)  by 
age  30  (AGEO),  with  an  average  diameter  of 
4.6  inches  (DBHO).  (The  data  deck  consisting 
of  32  cards  is  shown  in  fig.  2.) 


COLUMN  NUMBERS 


111  Ullll  122222222223333333333**********5 
1 2 345678901 2 3*5678901 2 34567890123*5678901 23*567890 


CARD  TYPE 
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2 

3 
* 
5 

3 
* 
5 

3 
* 
5 

3 
* 

5 

3 
* 
5 

3 
* 
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6 

*00 

30 

30 
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30 

30 
120 

30 

30 
120 

30 
30 
90 

30 
30 
90 

30 
30 
90 


10 
2000 


20 
*5 
80 

20 
45 
80 

20 
*5 
80 

20 
45 


20 
45 
80 

20 
45 
80 
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150 


850 
150 


850 
150 


850 
120 


350 
120 


850 
120 


80 
150 


80 
150 


80 
150 


Figure  2.  —  Data  deck  for  test  problem. 


An  Application  of  SPRYLD 

The  problem  described  below 
demonstrates  the  computations  made  by 


Yield  tables  produced  by  SPRYLD,  a  few 
of  which  are  reproduced  in  Appendix  2,  can 
assist  in  decision-making  in  many  ways. 
Money  yields  and  rates  earned  can  be  com- 
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Order  and  Contents  of  the  Data  Deck  for  Program  SPRYLD 


Card  Number  Variable 
type     of  cards  name 


Columns    Format    Description  of  variable 


1  NTSTS 
GIDE 

RINT 

1  COMCU 

COMBF 


3  1  per  JCYCL 

test 

MIX 


DSTINC 


1  per 
test 


AGEO 
DBHO 
DENO 
DSTY 


1  per 
test 


SITE 

THIN 

REGN(l)1 

VLLV(l) 


1-4  14  Number  of  tests  or  sets  of  yield  tables  to 

be  produced  (greater  than  or  equal  to  1). 

5-8  F4.0  Base  level  of  set  of  growing  stock  levels 

(equal  to  80.0  for  the  listing  in  Appen- 
dix 1). 

9-12  F4.0  Number  of  years  for  which  a  growth 

equation  makes  one  projection  (equal  to 
10.0  for  the  listing  in  Appendix  1). 

1-8  F8.3  Minimum  cut  in  merchantable  cubic 

feet  per  acre  to  be  included  in  total  yields 
(greater  than  or  equal  to  1.0). 

9  16  F8.3  Minimum  cut  in  board  feet  per  acre  to  be 

included  in  total  yields  (greater  than  or 
equal  to  1.0). 

1-4  14  Interval  between  intermediate  cuts.  A 

multiple  of  RINT. 

5-8  14  Number  of  stocking  levels  for  interme- 

diate cuts  to  be  examined  in  one  test. 
Equivalent  to  number  of  yield  tables 
produced  per  test  (greater  than  or  equal 
to  1). 

9-16  F8.3  Amount  growing  stock  level  (for  inter- 

mediate cuts  after  initial  thinning)  will 
be  increased  over  level  of  previous  yield 
table  in  a  test  if  MIX  is  greater  than  1. 
Leave  blank  if  MIX  equals  1. 
1-8  F8.3  Initial  age  in  years  to  be  shown  in  a 

yield   table.   Stand   age   when  first 
thinning  occurs  (greater  than  0.0). 
9-16  F8  3  Average  stand  d.b.h.  in  inches  just  prior 

to  initial  thinning  at  stand  age  AGEO 
(greater  than  0.0). 

17-24  F8.3  Number  of  trees  per  acre  just  prior  to 

initial  thinning  at  stand  age  AGEO 
(greater  than  0.0). 

25-32  F8.3  Lowest  growing  stock  level  for  inter- 

mediate cuts  after  initial  thinning  in  a 
test.  Level  will  increase  by  DSTINC  on 
the  second  and  subsequent  yield  tables 
in  a  test  if  MIX  is  greater  than  1.  Value 
of  DSTY  must  be  greater  than  0.0. 

33-40  F8.3  Site  index  for  the  stand  (greater  than 

0.0). 

41-48  F8,3  Growing  stock  level  for  initial  thinning 

at  age  AGEO  (greater  than  0.0). 

1-8  F8.3  Stand  age  at  which  first  regeneration 

cut  will  occur.  Must  be  greater  than  0.0 
as  REGN(l)  is  rotation  age  for  clear- 
cutting. 

9-16  F8.3  Percentage,  as  a  decimal,  of  growing 

stock  level  for  intermediate  cuts  to  be 
left  at  age  REGN(l).  Leave  this  and 
next  3  variables  blank  for  clearcutting. 
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Order  and  Contents  of  the  Data  Deck  for  Program  SPRYLD  (cont'd) 


Card     Number  Variable 

type      of  cards    name         Columns      Format    Description  of  variable 


REGN(2)1 

17-24 

F8.3 

Stand  age  at  which  second  regenera- 
tion cut,  if  any,  will  occur.  Final  cut  of 
2-cut  shelterwood  or  second  cut  of  3-cut 
shelterwood. 

VLLV(2) 

25-32 

F8.3 

Percentage,  as  a  decimal,  of  growing 
stock  level  left  after  first  regeneration 
cut  to  be  left  at  age  REGN(2).  Leave  this 
and  next  variable  blank  for  2-cut  shel- 
terwood. 

REGN(3)r 

33-40 

F8.3 

Stand  age  at  which  third  regeneration 
cut,  if  any,  will  occur.  Final  cut  of  3-cut 
shelterwood. 

1Values  for  ages  for  regeneration  cuts  must  equal  the  value  of  AGEO  plus  a  multiple  of  the 
value  ofRINT. 


puted  by  applying  thinning  costs  and  stum- 
page  values  to  the  volumes  given  in  the 
tables.  Stand  ages  at  culmination  of  mean 
annual  increment,  and  rates  earned  assist  in 
the  selection  of  rotations. 

For  the  situation  described  above,  yields 
and  numbers  of  precommercial  thinnings  are 
of  greatest  immediate  interest.  These  items 
are  summarized  in  tables  3  and  4  for  the  18 
yield  tables  produced.  Combinations  of  high 
initial  and  intermediate  subsequent  growing 
stock  levels  produce  the  greatest  volumes 
with  one  precommercial  thinning.  Additional 
comparisons  should  be  made  to  include  such 

Table  3. --Yields  in  board  feet,  including  com- 
mercial thinnings,  of  the  18  combinations  of 
initial  and  subsequent  growing  stock  levels1 

Initial  Subsequent  growing  stock  level 

growing  stock   

level  80  100  120 

  mbd  ft   

Two-cut  Shelterwood 

80  40.9  46.3  53.  4 

100  41.8  47.5  52.8 

120  42.1  47-9  52.2 

Three-cut  Shelterwood 

80  37.2  42.2  48.2 

100  38.4  43.6  48.2 

120  39.0  44.5  48.0 


factors  as  probable  thinning  costs,  cubic 
yields  from  thinnings  not  commercial  for 
board  feet,  and  the  average  size  of  tree 
produced.  As  expected,  the  current  crop 
produces  more  board  feet  in  150  years  if  cut  by 
two-cut  shelterwood  than  if  by  three-cut 
shelterwood.  The  latter  treatment  may, 
however,  get  the  next  crop  off  to  an  earlier 
start. 

Modifications  of  SPRYLD 

SPRYLD  can  be  modified  to  study  actual 
stands,  especially  to  determine  if  treatment  in 
unthinned  stands  is  justified,  and  to  add 
other  measures  and  variability  as  described 
for  PONYLD  by  Myers  (1971). 


Table  4. --Number  of  precommercial  thinnings  if 
each  of  the  18  combinations  of  initial  and 
subsequent  growing  stock  levels1  is  established 
as  specified  by  the  data  deck.     (Both  types 
of  cutting  gave  the  same  results) 


Initial  Subsequent  growing  stock  level 

growing  stock   —  

level  80  100  120 

80  12  2 

100  112 

120  11  2 


^ee  text  p.  3  for  description  of  growing  stock  ^ee  text  p. 3  for  descr i pt ion  of  growing  stock 
levels.  levels. 
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Description  of 
Program  SPRCHK 

Program  SPRCHK  —  used  to  calculate 
volumes,  volume  conversion  factors,  and 
diameters  and  heights  for  different  com- 
binations of  stand  variables  not  available 
from  growth  prediction  plots  —  is  the  same  as 
Myers  (1971)  Program  PONCHK,  except  that 
the  equations  used  to  compute  volumes  are  for 
spruce  (Myers  and  Edminster  1972). 

Modifications  Needed  to  Use 
Spruce-Fir  in  Myers'  (1974) 
Program  TEVAP2 

Subroutine  WORKGP  in  Myers  (1974) 
TE  VAP2  permits  the  program  to  be  used  with 
other  species.  To  use  spruce-fir,  replace  the 
dummy  continue  statement  number  4  with  a 
CALL  statement  that  will  call  subroutine 
ESSF,  the  species-specific  statements  for 
spruce-fir  to  be  used  with  TEVAP2.  SPNUM 
(I)  equal  to  4  will  then  designate  subroutine 
ESSF.  The  subroutine  organization  is  shown 
in  Appendix  3. 

Basic  information  used  in  ESSF  has  been 
described  earlier,  with  the  exception  of  one 
equation  (ADD)  which  is  included  to  estimate 
merchantable  cubic  volume  obtained  as  a 
byproduct  of  saw-log  cuts.  ADD  is  computed 
from  stand  diameter  as  follows: 

ADD  =  1.54375  +  1L43324 

D 
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S     =  0.1168,   R2  =  0.8517 


where        ADD  =  cubic  foot  volume  of  saw 
logs  (hundreds  of  cubic 
feet  per  m  board  feet). 
D  =  average  stand  diameter. 
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Appendix  I:  Listing  of  Program  SPRYLD 


PROGRAM  SPRYLD 
1 ( INPUT, OUTPUT, TAPE5*INPUT.TAPE6=0UTPUT) 

TO  COMPUTE   AND  PRINT   YIELD  TABLES  FOR  MANAGED  EVEN-AGED   STANDS  OF 
ENGELMANN   SPRUCE   AND   SUBALPINE  FIR. 

DEFINITIONS   OF  VARIABLES. 

ADDHT   =    INCREASE   OR  DECREASE   IN   AVERAGE   STAND  HEIGHT   BY  THINNING. 

AGEO  -   INITIAL  AGE   IN  YIELD  TABLE. 

BASC   =   BASAL   AREA  CUT   PER  ACRE. 

BASO  =   BASAL   AREA   PER   ACRE   BEFORE  THINNING. 

BAST  =   BASAL   AREA  PER   ACRE  AFTER  THINNING. 

BDFC   =   BOARD  FEET  CUT   PER  ACRE. 

BDFO  =   BOARD  FEET   PER   ACRE   BEFORE  THINNING. 

BDFT  =   BOARD  FEET  PER   ACRE   AFTER  THINNING. 

CFMC  =   MERCHANTABLE  CU.   FT.   CUT   PER  ACRE. 

CFMO  -   MERCH.   CU.   FT.    PER  ACRE   BEFORE  THINNING. 

CFMT   =   MERCH.   CU.   FT.    PER   ACRE  AFTER  THINNING. 

COMBF  =   MINIMUM  COMMERCIAL  CUT,    BOARD  FEET. 

COMCU  =   MINIMUM  COMMERCIAL   CUT,   CU.  FT. 

OBHO  =  AVERAGE   STAND  D.B.H.    BEFORE  THINNING. 

DBHT  =  AVERAGE   STAND   D.B.H.   AFTER  THINNING. 

DENC  =  TREES  CUT  PER  ACRE, 

DENO  =   TREES   PER   ACRE   BEFORE  THINNING. 

DENT   =   TREES   PER   ACRE   AFTER  THINNING. 

DIED  =   PERCENTAGE,    AS   A  DECIMAL,   OF  TREES   THAT  DIE   DURING  PERIOD 
RINT. 

OLEV  =  GROWING  STOCK   LEVEL   FOR   INTERMEDIATE   CUTS   AFTER  FIRST. 
DSTINC   =   AMOUNT   GROWING  STOCK  LEVEL   FOR    INTERMEDIATE   CUTS   WILL  BE 

INCREASED  OVER   LEVEL   OF   PREVIOUS   YIELD  TABLE    IN  A  TEST. 
D STY  =  LOWEST  VALUE  OF  DLEV  USED   IN  A  TEST. 

GIDE   =  BASE   FOR  GROWING  STOCK   LEVELS,    80.0   IN   EXAMPLE  SHOWN. 

HTSO  =   TREE  HEIGHT  BEFORE  THINNING. 

HTST  =   TREE  HEIGHT  AFTER  THINNING. 

JCYCL  =   INTERVAL   BETWEEN   INTERMEDIATE  CUTS. 

JSBD  =   SUM  OF   BOARD  FEET   FROM  ALL   CUTS  WITH  YIELD  OF  COMBF  OR 
LARGER. 

JSMC   =   SUM  OF   MERCH.   CU.   FT.   FROM   ALL   CUTS  WITH  YIELD  OF  COMCU  OR 
LARGER. 

JSTF   =   SUM  OF   TOTAL   CU.   FT.   FROM   ALL  CUTS. 

MIX  =   NUMBER   OF   STOCKING  LEVELS   EXAMINED   PER  TEST. 

NTSTS  =  NUMBER  OF   TESTS  PER  BATCH. 

PRET  =   PERCENTAGE   OF   TREES   RETAINED  AFTER  THINNING. 
REGN(I)    =  STAND   AGE  WHEN  REGENERATION   CUT   I  OCCURS. 
RINT  =  NUMBER  OF   YEARS  FOR  WHICH   PROJECTION   IS  MADE. 
ROTA  =  FINAL   AGE    IN  YIELD  TABLE. 
SITE   =   SITE  INDEX. 

THIN  =   GROWING  STOCK   LEVEL   FOR   INITIAL  THINNING. 
TOTC  =   TOTAL   CUBIC   FEET  CUT   PER  ACRE. 
TOTO  =   TOTAL  CUBIC   FEET  PER   ACRE   BEFORE  THINNING. 
TOTT  =   TOTAL  CUBIC   FEET   PER   ACRE   AFTER  THINNING. 

VLLV(I)    =   PERCENT   OF   PREVIOUS  DLEV  TO  BE   LEFT   AT  REGN(I),  ENTERED 
AS   A  DECIMAL. 


COMMON  BA, BAST, CUFT.DBHO.OBHT, DENO, FCTR.H IT E, GIDE, PRET, PR QD.PE ST, S 
1 TAND , VDM 
01  MENS  I  ON  REGN(3),VAR(10),VLLV(2) 

DO   1  J=l,10 
1   VAR(J)    =  0.0 

READ  NUMBER   OF   TESTS,   BASE   OF   GROWING   STOCK  LEVELS,    AND  LENGTH  OF 
PROJECTION  PERIOD   FROM  CARD  TYPE  ONE. 

READ   (5,5)  NTSTS, GIDE, RINT 
5   FORMAT  (I4.2F4.0) 

IFINTSTS   .LE.   0)   GO  TO  170 
IF (GIDE   .LE.   0.0)   GO  TO  170 
VARI5)    =  RINT 

READ   MINIMUM  COMMERCIAL  CUTS   FOR   COMPUTATION   OF   COLUMN   TOTALS  FROM 
CARD   TYPE  TWO. 

READ    (5,10)   COMCU, COMBF 
10  FORMAT  (10F8.3) 
VARI8)    =  COMCU 
VARI9)    =  COMBF 

EXECUTE   PROGRAM  ONCE  FOR   EACH  SET   OF    INITIAL   VALUES  OF  INTEREST. 

DO   160   1=1, NTSTS 
JTEM  =  0 

READ  CUTTING   INTERVAL   AND  LEVELS   PER  TEST  FROM  CARD  TYPE  THREE. 

READ    (5,15)    JCYCL, MIX, OSTINC 
15   FORMAT    1214, F8. 3) 

IFIJCYCL    .LE.   0   .OR.    MIX   .LE.   0)    GO  TO  170 
JTEM   =  JCYCL 

READ   INITIAL   STANO  VALUES  FROM  CARD   TYPE   FOUR • 


READ   (5,10)    AGEO, DBHO , DENO, DSTY, SITE, THIN 


VAR ( 1  ) 

=  AGEO 

VARI2) 

=  DBHO 

VAR(3) 

=  DENO 

VAR (4) 

=  DSTY 

VAR (6 ) 

=  SITE 

VAR ( 7 ) 

=  THIN 

C   READ   SILVICULTURAL  CONTROLS  FROM  CARO  TYPE  FIVE. 

C 

READ    (5,10)  REGN(l),VLLV(l),REGN(2),VLLV(2),REGN(3) 
VAR(IO)   =  REGN(l) 
DO  20  L=l,10 

IF(VAR(L)    .LE.   0.0)   GO  TO  170 

20  CONTINUE 
DLEV  =  0.0 
DO  25  NA*1,3 

L  =   4  -  NA 

I F ( REGN ( L )    .EO.   0.0)    GO  TO  25 
ROTA   =  REGN(L) 
GO  TO  30 
25  CONTINUE 


PROVIDE   FOR   SEVERAL   GROWING  STOCK  LEVELS   PER  TEST. 


30  DO   160  H=1,MIX 
ADDHT  =  0.0 
BDFO  =  0.0 
BOFT  =  0.0 
CFMO   =  0.0 
CFMT  =  0.0 
HTCUM  =0.0 
JSBD  =  0 
J  SMC  =  0 
JSTF   =  0 
TEM  =  M 
DLEV  =  DSTY 
BASO  =  DENO 


(TEM  -  1.0) 
0.0054542  * 


*  OSTINC 
DBHO  »  DBHO 


OBTAIN  AVERAGE  HEIGHT   AND  VOLUMES   PER  ACRE. 

 STATEMENTS  FOR  HTSO  AND   IF  STATEMENT   ARE  SPECIES-SPECIFIC. 


IF ( AGEO   .GE.    100.0)   GO  TO  35 
HTSO  =  -13.71751  ♦  0.15087   *  SITE 
1   0.01371   *  AGEO  *  SITE  -  0.00006  * 
GO  TO  40 

35  HTSO  =  0.91859  -   100.43601   /  AGEO 
1   40.08154  *  ALOGIO(SITE)    /  AGEO 

HTSO  =   10.0   **  HTSO 
40  HITE   =  HTSO 

BA  =  BASO 

STAND  =  DENO 

VDM   =  DBHO 

CALL  SPRVOL 

TOTO  =  CUFT 

BDFO  =  CUFT  »  PROO 

CFMO  =   CUFT   *  FCTR 

REST  =  THIN 


ENTER   LOOP  FOR  REMAINING  COMPUTATIONS   AND  PRINTOUT. 

DO  130  K=l,100 

CHANGE   STANDARDS    IF   A  REGENERATION  CUT    IS   DUE • 

43   I F ( AGEO   .GE.   ROTA)   GO  TO  60 

IFIAGEO   .LT.    REGN ( 1 ) )   GO  TO  55 
IF ( AGE 0   . NE.    REGN(ll)    GO  TO  45 
DLEV  =  DLEV  *  VLLV(l) 
REST  =  OLEV 

JCYCL  =   REGN ( 2 )    -  REGN ( 1 )    *  0.5 

GO   TO  55 
45  DLEV  =  DLEV  *  VLLVI2) 
REST  =  DLEV 
JCYCL   =  REGN ( 3  > 


»  0.00126  »  AGEO  *  AGEO 
AGEO   *  AGEO  *  SITE 


0.62318  «  ALOGIO(SITE)  ♦ 


REGNI2)   +  0.5 


INCREASE   D.B.H.   BY  THINNING  AND  COMPUTE  POST-THINNING  VALUES. 


55  CALL  SPRCUT 

I F ( PRE T  .GE.  100.0)  GO  TO  56 
JOENT  =  (BAST  /  (0.0054542  * 
DENT  =  JDENT 

BAST  =  0.0054542  *  DBHT  *  DBHT 
IF ( BAST   .LT.   BASO)    GO  TO  58 

56  BAST   =  BASO 
HTST  =  HTSO 
DENT  =  DENO 
JDENT   =   DENO  +  0.5 
DBHT   =  DBHO 

TOTT  =  TOTO 
BDFT  =  BDFO 
CFMT  =  CFMO 
GO  TO  60 


DBHT   *   DBHT ) )   ♦  0.5 
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SORT [ 1.0   /   PR  6  T  1 


£  STATEMENT  FOR  AOOHT   IS  SPECIES-SPECIFIC. 

C    S3  ADDHT   =  -   1.81392   ♦  17.77922 
HTCUM  -  HTCUM  ♦  ADOHT 
HTST   =  HTSO  ♦  ADOHT 
STAND  =  DENT 
VDM  =  08HT 
HI  TE   =  HTST 
8A  =  RAST 
CALL  SPRVOL 
TOTT  =  CUFT 
BDFT  =  CUFT  *  PROD 
CFMT  =  CUFT   »  FCTR 

C 

I  CHANGE  MODE   AND  ROUND  OFF   FOR  PRINTING. 


60  JAGEO 
JSITE 
JDENO 
JHTSO 
JTOTO 
JTOTO 
J8AS0 
JCFMO 
JCFMO 
JBOFO 
JBDFO 
JHTST 
JTOTT 
JTOTT 
JCFMT 
JCFMT 


AGED 
SITE 
DENO  ♦ 
HTSO  ♦ 
(TOTO 
JTOTO 
BASO  + 
(CFMO 
JCFMO 
(  8DF0 
JBDFO 
HTST  ♦ 
I  TOTT 
JTOTT 
(CFHT 
JCFMT 


0.5 
0.5 

'  0.1)  < 
'  10 
0.5 

'   0.1)  < 
'  10 
'  0.01) 
•  100 
0.5 

■0.1)  t 
'  10 

'  0.1)  < 
<  10 


IF ( JCFMT  . GT . 
JflDFT  =  [BDFT 
JBDFT  =  JBDFT 
IFIJ8DFT  .GT. 
J8AST  =  BAST 
JOENC  =  JDENO  -  JDENT 
J8ASC  =  JBASO  -  J8AST 
JTOTC  =  JTOTO  -  JTOTT 
JCFMC  =  JCFMO  -  JCFMT 
LE 


JCFMO)    JCFMO  =  JCFMT 

*  0.01)  +0.5 

*  100 

JBDFO)  JBOFO  =  JBDFT 
0.5 


IF( JCFMC 


0)    JCFMC  =  0 


JBDFC   =   JBDFO  -  JBDFT 

IFUBOFC    .LE.  0)    JBDFC  =  0 

C 

;  SUM  PERIODIC   CUTS  FOR  LAST  LINE   OF   YIELD  TABLE. 

c 

IFIAGEO   .GE.    ROTA)    GO  TO  70 
JSTF  =   JSTF   ♦  JTOTC 
CFMC  =  JCFMC 

IFICFMC   • LT.   COMCUI    GO  TO  65 
JSMC  =  JSMC  +  JCFMC 
65  BDFC  =  JBDFC 

IFIBDFC   .LT.   COMBF )   GO  TO  70 
JSBO  =   JSBD  ♦  J80FC 

C 

C  WRITE  HEADINGS   FOR  YIELD  TABLE. 

C 

70  IFIK  .GE.   2)   GO  TO  92 


 CHANGE   TABLE   HEADING   FOR  OTHER  SPECIES. 

WRITE    16,80)    JSITE, JCYCL 
80  FORMAT    ( 1H1 , / / , 27X , 82 HY I  ELDS   PER   ACRE  OF   MANAGED,    EVEN-AGED  STANDS 
1   OF   ENGELMANN   SPRUCE   AND  SUBALPINE   FIR/1HC48X,  11HSITE    INDEX  ,13,1 
2H,,I4,19H-YEAR  CUTTING  CYCLE) 
WRITE    16,82)    THI N, DLEV 
82  FORMAT    ( 1H0.41X.26HTHINNING  LEVELS=   INITIAL   -,F6.0,14H,  SUBSEQUENT 
1  -,F6.0) 
WRITE    16, 8*1 

84  FORMAT    (1H0,25X,38HENTIRE  STAND  BEFORE   AND  AFTER   T  HINN  ING  ,  28x  ,  26eiP 
1ERI0DIC   INTERMEDIATE  CUTS) 
WRITE  (6,86) 

86  FORMAT    (1HO,9X,5HSTAND,10X,5HBASAL,3X,7HAVERAGE,2X, 7HAVERAGE , 3X.5H 
1  TOTAL , 3X .9HMERCHANT-, 3X, 9HS AWT  I MBER , 9X , 5H BASAL ,  4X, 5HT0TAL ,3X,9HMER 
2CHANT-, 3X , 9HSAWTIMBER ) 
WRITE  (6,88) 

88  FORMAT   (1H   ,  1  OX  ,  3HAGE  ,  4X  ,  5HTR  EES  ,  3X ,  4HAR  F  A ,  4X  ,  6HD .  B  .H  .  ■  3X  ,  6HHE  I  GHT 
1  ,2X,6HV0LUME, 2X, 11HABLE  V0LUM£,4X»6HV0LUME»3X,5HTREES»3X,4HAREA,3X 
2.6HV0LUME ,2X, 11HABLE   VOLUME , 4X , 6HV0L UME ) 
WRITE  (6,90) 

90  FORMAT    (1H    , 8X , 7H ( YEARS ) , 3X , 3HN0 . , 3X , 6HSQ . FT . , 4X , 3H I N . , 6X , 3HF T. , 4X 


1 .6HCU.FT. , 5X, 6HCU.FT. ,6X, 6HB0.FT. ,4X, 3HN0. 
2,5X,6HCU.FT. ,6X,6HBD.FT. ) 


3X.6HSQ.FT. ,2X,6HCU.FT. 


C 

C  WRITE   TABLE   ENTRIES   OF   DIAMETER,    VOLUMES ,  ETC. 

92  WRITE    (6,94)    JAGEO, JDENO, JBASO, OBHO, JHTSO, JTOTO, JCFMO, JBDFO 
94  FORMAT  (1H0,9X,I4,4X,I5,2X,I4,5X,F5.1,5X,13,4X,I5,6X,I5,6X,I6) 
IFIAGEO   .GE.   ROTA)    GO   TO  135 

WRI TE    (6,96)    JAGEO, JDENT, JB AST, DBHT, JHTST, JTOTT, JCFMT, JBDFT, JDENC, 
1 JBASC , JTOTC, JCFMC, J80FC 
96  FORMAT    (  1H   ,  9X  ,  I  4 ,  4X ,  I  5  ,  2  X ,  I  4 ,  5X  ,  F5  .  1 ,  5X  ,  I  3  ,  4X  ,  I  5  ,  6X  ,  I  5  ,  6  X  ,  I  6  ,4X  ,  I 
15,3X,I3,5X,I4,6X,I4,8X,I5I 

C  COMPUTE   VALUES  FOR  EACH  PERIOD.   THIN   AS  SPECIFIED. 

IRINT  =  R I  NT 

IK  =   JCYCL  /  IRINT 

DO  120  L=1,IK 

AGEO  =   AGEO  +  RINT 

IFIAGEO   .GT.   ROTA)   GO  TO  135 

C 

C  COMPUTE   NEW  D.B.H.   BEFORE   THINNING  AND  ROUND  OFF   TO   0.1  INCH. 


--STATEMENT   FOR   DBHO   IS  SPECIES-SPECIFIC. 

DBHO   =   1.62917  ♦   1.03371   *   DBHT   ♦  0.01304  •   SITE  -   0.90669  • 
1  ALOGICIBAST) 
IOBHO  =   OBHO   »   10.0   *  0.5 
DBHO  =    I DBHO 
DBHO   =   DBHO  *  0.1 


-STATEMENT  FOR  OIEO  IS  SPECIES-SPECIFIC.  CHANGE  10.0  IN  IF 
-STATEMENT    IF   DIED   STATEMENT   APPLIES   TO  LARGER  TREES. 


I F ( DBHT   .GE.    10. 0)    GO  TO  100 

DIED  =    1-0.03967   +  0.00382   *   DBHT   »   BAST)   *  0.01 
IFIDIED   .LT.   0.0)    DIEO  =  0.0 
DENO   =   DENT   «    11.0  -  OIED) 
MNK  =  OENO  ♦  0.5 
DENO  =  MNK 
GO  TO  105 
100  DENO  =  DENT 

105  BASO  =  DENO  *  10.0054542  »  DBHO  *  DBHO) 
OBTAIN   AVERAGE  HEIGHT  AND  VOLUMES   PER  ACRE. 

 STATEMENTS  FOR  HTSO  AND   IF   STATEMENT   ARE  SPECIES-SPECIFIC. 


IFIAGEO   .GE.    100.0)    GO  TO  110 

HTSO   =  -13.71751   +  0.15087   *   SITE   ♦  0.00126  *  AGEO  «   AGEO  ♦ 
1   0.01371   *  AGEO  *  SITE  -  0.00006  *   AGEO  *  AGEO  *  SITE 
GO  TO  115 

110  HTSO   =   0.91859  -   100.43601   /   AGEO   ♦  0.62318   ♦   AL0G10ISITE)  ♦ 
1    40.08154   «   ALOG10ISITE)    /  AGEO 

HTSO  =  10.0  «*  HTSO 
115  HTSO  =   HTSO  +  HTCUM 

STAND  =  DENO 

VDM  =  DBHO 

H I TE   =  HTSO 

BA   =  BASO 

CALL  SPRVOL 

TO  TO   =  CUFT 

BDFO  =  CUFT   *  PROD 

CFMO  =   CUFT   *  FCTR 

TEST   IF   REGENERATION  CUT    IS  DUE. 


DO  118  KU=1,3 

IFIAGEO   .EO.   REGN(KU))   GO  TO  43 
118  CONTINUE 

CHANGE   MODE   AND  ROUND  OFF   FOR  PRINTING. 


IFIL   .EO.    IK)   GO  TO  125 

KOENO  =   DENO  +  0.5 

KAGEO  =  AGEO 

KHTSO  =  HTSO  +  0.5 

KBASO  =  BASO  +  0.5 

KTOTO  =    (TOTO  *  0.1)    +  0.5 

KT9T0  =   KTOTO  *  10 

KCFMU  =   I CFMO  *  0.11   +  0.5 

KCFMO  =   KCFMO  *  10 

KBDFO  =    ( BDFO  *  0.01)    ♦  0.5 

KBDFO  =   KBOFO  *  100 


WRITE   VALUES   FOR  THE   PERIOD   IF   THINNING   IS  NOT  DUE • 

WRITE    (6,94)    KAGEO, KOENO, KBASO, DBHO , KHTSO, KTOTO, KCFMO, KBDFO 

DBHT  =  DBHO 

BAST  =  BASO 

OENT  =  DENO 
120  CONTINUE 
125   REST  =  DLEV 
130  CONTINUE 

ADD  FINAL  CUTS   TO  TOTAL   YIELDS   AND  WRITE  TOTAL  YIELDS. 

135   JSTF   =   JSTF   ♦  JTOTO 
CFMO  =  JCFMO 

IFICFMO   .LT.   COMCU)    GO  TO  140 

JSMC  =  JSMC  +  JCFMO 
140   BDFO  =  JBDFO 

I F ( BDFO   .LT.   COMBF )    GO  TO  145 

JSBD  =  JSBD  *  JBDFO 
145   WRITE    (6,150)    JSTF, JSMC, JSBD 

150   FORMAT    I 1H0 , / , 67X , 1 2HT0T AL   Y I  EL  OS , 18X , I  6 , 4X , 1 6 , 6X , I  7 ) 
WRITE    (6,155)    COMCU, COMBF 

155  FORMAT    ( 1H0 , / / , 1 1 X ,44HM I N I  MUM  CUTS   FOR    INCLUSION    IN  TOTAL   YIELDS — 
1.F6.0.15H  CUBIC   FEET    AND,F7.0,11H   BOARD  FEET ) 

WRITE  16,156) 

156  FORMAT    ( 1H0 , 1  OX , 66HME RCH.   CU.   FT.   -  TREES   5.0    INCHES   D.B.H.   AND  LA 
1RGER   TO  4-INCH  TOP. 1 

WRITE  (6,157) 

157  FORMAT    ( 1H0, 10X , 59HBD.   FT.   -   TREES   8.0   INCHES   D.B.H.    AND  LARGER  TO 
1    6- INCH  TOP.) 

PREPARE   FOR  NEXT   TABLE   OF  THE  TEST. 

AGEO   =  VAR(l) 

DBHO  =  VAR(2) 

OENO  =   VAR (3 ) 

JCYCL  =  JTEM 
160  CONTINUE 

GO   TO  200 
170   WRITE  (6,175) 

175   FORMAT    (1H1,///,10X,66HEXECUTI0N   STOPPED   BECAUSE   OF  NEGATIVE   OR  ZE 

1R0   ITEM  ON  A  DATA  CARD.) 
200  CALL  EXIT 

END 
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SUBROUTINE  SPRCUT 

TO  ESTIMATE   INCREASE    IN  AVERAGE   D.B.H.   DUE  TO  THINNING. 

COMMON  BA, BAST, CUFT, DBHO, DBHT»DENO»FCTR,HITE»GIDE, PRET, PRCD, REST, S 
1 TANO  t  VDM 

IFIDBHO   .IT.   9.*)    GO  TO  30 

COMPUTE   O.B.H.    IF   OBHO   IS   LARGE   ENOUGH  FOR   BASAL   AREA   TO  REMAIN  CONSTANT. 

PRET   =  100.0 

DO  21   KJ=l,100  C 


-STATEMENTS  FOR  DBHE   AND  PDBHE   ARE  SPECIES-SPECIFIC. 


I F  ( PRET    .LT.    50.0)   GO   TO  5 
DBHE   =  0.02666  *   1.30655   *  DBHO 
GO  TO  11 

5   PDBHE   =   0.33206   ♦  0.983*6   *   AL0G10 ( DBHO ) 
10.0  *»  PDBHE 
=   DBHE     *   10.0   ♦  0.5 
I  DBHE 

DBHE   *  0.1 
DENO  *   PRET   *  0.01 
=   OENE   ♦  0.5 
NDENE 

0.0054542   *   DBHE   *  DBHE   *  DENE 
=   BASE   *   10.0  ♦  0.5 
NBASE 

BASE   *  0.1 

0. 00545*2   *  DBHE  *  DBHE 
TEM  -  BASE  -  REST 

IFIKJ   .£0.    1    .AND.   TEM   .LT.   0.0)    GO  TO  90 
I F I TEM   .LE.   TMPY)   GO   TO  70 
IFITEM   .LT.   4.0)    GO   TO  20 
PRET  =   PRET  -  1.0 
GO  TO  21 

20  PRET  =   PRET  -  0.3 

21  CONTINUE 
GO   TO  70 

COMPUTE   D.B.H.    IF   BASAL   AREA   INCREASES  WITH  D.B.H. 

30  PRET  =  40.0 

IFIDBHO   .GT.   7.0)    PRET  =  70.0 
DO  65  J=l,100 


0.00306   »   PRET   *  DBHO 

0.14710   *  ALOGIO(PRET) 


DBHE  = 
11    I DBHE 
DBHE  = 
DBHE  = 
DENE  = 
NDENE 
DENE  = 
BASE  = 
NBASE 
BASE  = 
BASE  = 
TMPY  = 


•STATEMENTS   FOR   DBHE   AND  PDBHE   ARE  SPECIES-SPECIFIC. 


IF!  PRET   .GE.    50.0)   GO   TO  40 

PDBHE   =   0.33206  +   0.98346   *   AL0G10 ( DBHO )    -  0.14710  *   ALOG 1 C ( PRE T ) 
DBHE   =   10. 0   **  PDBHE 
GO  TO  45 

40  DBHE   =   0.02666  ♦   1.30655   *  DBHO  -  0.00306   *   PRET   *  OBHO 
45   I DBHE   =   DBHE   *   10.0  +  0.5 

DBHE   =   I DBHE 

DBHE   =   DBHE   «  0.1 

OENE   =  DENO  *   (PRET   *  0.011 

NDENE   '   OENE   +  0.5 

DENE   =  NDENE 

BASE   =  0.0054542   *   DBHE   *   DBHE   *  DENE 
NBASE   =   BASE   *   10.0   +  0.5 
BASE   =  NBASE 
BASE   =   BASE   «  0.1 


-CHANGE  STATEMENTS  FOR  BREAK,  BUST  AND  FIRST  3  STATEMENTS  FOR  DBHP 
-IF   OTHER  GROWING   STOCK   LEVEL   BASE   THAN   TABLE   1    IS  USED. 


BREAK  =   49.9   *  REST   /  GIDE 
I F ( BASE   .GT.   BREAK)   GO  TO  50 

DBHP   =    (GIDE   /   REST)    *   (0.08682   *   BASE)    +  0.94636 
GO  TO  52 

50  BUST   =   66.2  »    (REST   /  GIDE) 
I F ( BASE    .GT.    BUST )   GO   TO  51 

DBHP   =    (GIDE   /   REST)    *   (0.10938   *   BASE)   -  0.17858 
GO  TO  52 

51  TMPY  =   BASE   *   (GIDE   /  REST) 
TEM  =   TMPY  *  TMPY 

DBHP   =   19.04740   *   TMPY  -  0.26673   *  TEM   +   0.0012539   *  TEM 
1   -  448.76833 
I F ( TMPY   .GT.   GIOE)    DBHP   =   DBHO  +  0.8 

52  I DBHP  =   DBHP   *   10.0   *  0.5 
DBHP  =   I DBHP 

DBHP  =   DBHP   *  0.1 
IFIOBHP  -   DBHE)  60,70.61 

60  PRET  =   PRET   *  1.02 

I F ( PRET  .GT.   100.0)   GO  TO  90 
GO  TO  65 

61  PRET  =  PRET  *  0.98 
65  CONTINUE 

70  DBHT   =  DBHE 

COMPUTE   POST-THINNING   BASAL  AREA. 


-CHANGE  TWO  IF  STATEMENTS  AND  STATEMENTS  FOR  SOFT  IF  DIFFERENT 
-GROWING   STOCK  LEVEL   BASE    IS  USED. 


77   BAST   =  REST 
80  RETURN 
90   PRET  =  100.0 
RETURN 

END 


SUBROUTINE  SPRVOL 

TO  COMPUTE   VOLUMES   PER  ACRE   IN  VARIOUS  UNITS. 

COMMON  BA,8AST,CUFT, DBHO, DBHT ,DENO,FCTR,HITE,GI0E,PRET,PR00,REST,S 
IT AND, VDM 

FCTR  =  0.0 
PROD  =  0.0 

COMPUTE   TOTAL   CUBIC   FEET   PER  ACRE. 

D2H   =   VDM   *  VOM  *  HITE 

 STATEMENTS  FOR  CUFT   AND   IF   STATEMENT   ARE  SPECIES-SPECIFIC. 


IF1D2H   .GE.    22500.0)    GO  TO  5 

CUFT  =    (0.0023575   *  02H  -   0.00524   *   BA   ♦   0.34480)    *  STAND 
GO  TO  10 

5  CUFT  =    (0.0020103   *  D2H  -  0.01052   *   BA   +  7.60196)   *  STAND 
10    IF (VDM   .LT.    5.0)    GO  TO  40 

OBTAIN  CONVERSION  FACTORS  FOR  MERCH.  CU.  FT.  -  VOLUMES  TO  4.0-INCH  TOP 
IN   TREES   5.0    INCHES   D.B.H.   AND  LARGER. 


-STATEMENT   FOR   FCTR   AND   IF   STATEMENT   ARE  SPECIES-SPECIFIC. 


FCTR   =   0.82375   ♦   3.45569  /   VDM  +   0.00013   *   VDM  *   VDM  -  28.16783  / 
1    ( VDM  *   VDM ) 
I F ( FCTR   .GT.    0.99)    FCTR  =  0.99 
25    IFIVDM   .LT.   8.0)    GO  TO  40 


OBTAIN  CONVERSION  FACTORS  FOR  BD.  FT. 
TREES   8.0   INCHES  D.B.H.    AND  LARGER. 


VOLUMES  TO   6.0-INCH  TOP  IN 


-STATEMENTS  FOR   PROD  AND   IF   STATEMENT   ARE  SPECIES-SPECIFIC. 


IFIVDM   .GE.    16.5)   GO  TO  30 

PROD  =  4.59159  -  214.06370  /  (VDM  *  VDM)  t  0.39782  *  ALOGiCIBAI 
GO  TO  40 

30  PROD   =   8.59422  -   19.54507   /   SORT(VOM)    +  0.44432   *  ALOGICIBAI 
40  RETURN 
END 


IF  I DBHT   .GT.    5.0)    GO   TO  75 

SOFT  =   11.58495   *   DBHT   -  11.09724 

GO   TO  76 

75  I F I DBHT   .GE.    10.0)    GO   TO  77 
TEM  =  DBHT   »  DBHT 

SOFT   =   7.76226   *  DBHT   *0. 85289  *   TEM  -0.07952   *   TEM  *   DBHT-3. 45624 

76  BAST  =    (REST   /  GIDE)    *  SOFT 
GO  TO  80 


MA- 


Appendix  II:  Output  of  SPRYLD 

Two-Cut  Shelterwood 

YIELDS   PER  ACRE   OF  MANAGED,    EVEN-AGED  STANDS   OF   ENGELMANN   SPRUCE   AND   SUBALP I NE  FIR 
SITE    INDEX     80,      30-YEAR  CUTTING  CYCLE 
THINNING  LEVELS*    INITIAL   -       80.,    SUBSEQUENT   -  80. 
ENTIRE   STAND  BEFORE   AND  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS  > 

TREES 
NO. 

BASAL 
AREA 
SO. FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT . 

SAWT I MBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

30 
30 

850 
314 

94 

52 

4.5 
5.5 

28 

29 

1010 
680 

340 
340 

0 
0 

536 

42 

330 

0 

0 

40 

311 

78 

6.8 

38 

1260 

900 

0 

SO 

305 

106 

8.0 

45 

2030 

1650 

4200 

60 
60 

295 
149 

133 
80 

9.1 
9.9 

53 
53 

2920 
1810 

2530 
1610 

8300 
5700 

146 

53 

1110 

920 

2600 

70 

145 

99 

11.2 

59 

2510 

2310 

9200 

80 

145 

122 

12.4 

65 

3370 

3150 

13600 

90 
90 

145 

55 

146 

80 

13.6 
16.3 

70 
71 

4350 
2440 

4120 
2340 

18700 
11100 

90 

66 

1910 

1780 

7600 

100 

55 

95 

17.8 

75 

3000 

2910 

14500 

110 

55 

112 

19.3 

79 

3610 

3520 

18300 

120 

55 

130 

62 

20.8 
25.2 

82 
84 

4290 
2050 

4200 
20  30 

22500 
1 1  300 

37 

68 

2240 

2170 

130 

18 

72 

27.1 

87 

2420 

2400 

13700 

140 

18 

83 

29.0 

89 

2830 

2810 

16500 

150 

18 

94 

30.9 

91 

3280 

3250 

TOTAL 

19500 
YIELOS 

8870 

8120 

40900 

MINIMUM 

CUTS 

FOR   INCLUSION  IN 

TOTAL  YIELDS- 

400.  CUBIC  FEET 

AND  2000. 

BOARD 

FEET 

MERCH. 

CU.  FT 

.   -  TREES 

5.0  INCHES  D.B.H. 

AND  LARGER  TO  4-INCH 

TOP. 

BD.   FT.   -   TREES   8.0   INCHES   D.B.H.   AND  LARGER  TO  6- INCH  TOP. 

YIELDS   PER  ACRE  OF  MANAGED,    EVEN-AGED  STANDS  OF   ENGELMANN   SPRUCE   AND   SUBALPINE  FIR 
SITE   INDEX     80,     30-YEAR  CUTTING  CYCLE 
THINNING  LEVELS*    INITIAL   -       80.,    SUBSEQUENT  -  120. 
ENTIRE   STAND   BEFORE   AND  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

30 
30 

850 
314 

94 
52 

4.5 
5.5 

28 
29 

1010 

680 

340 
340 

0 
0 

536 

42 

330 

0 

0 

40 

311 

78 

6.8 

38 

1260 

900 

0 

50 

305 

106 

8.0 

45 

2030 

1650 

4200 

60 
60 

295 
228 

133 
119 

9.1 
9.8 

53 
53 

2920 
2660 

2530 
2360 

8500 
8  500 

67 

14 

260 

170 

0 

70 

218 

141 

10.9 

59 

3520 

3220 

12800 

80 

218 

171 

12.0 

65 

4670 

4350 

18700 

90 
90 

218 

92 

204 
121 

13.1 
15.5 

70 
71 

5990 
3650 

5640 
3500 

25500 
16500 

126 

83 

2340 

2140 

9000 

100 

92 

142 

16.8 

75 

4570 

4400 

21800 

110 

92 

164 

18.1 

79 

5320 

5160 

26500 

120 
120 

92 
32 

189 
96 

19.4 
23.4 

82 
83 

6240 
3150 

6080 
3120 

32300 
17100 

60 

93 

3090 

2960 

15200 

130 

32 

110 

25.1 

86 

3700 

3670 

20700 

140 

32 

125 

26.8 

89 

4310 

4270 

24800 

150 

32 

142 

28.5 

91 

4960 

4910 

29200 

TOTAL 

YIELDS 

10980 

10010 

53400 

MINIMUM 

CUTS 

FOR   INCLUSION  IN 

TOTAL   YIELDS —     400.  CUBIC 

FEET 

AND  2000. 

BOARD 

FEET 

MERCH. 

CU.  FT 

.   -   TREES   5.0   INCHES   D.B.H.   AND  LARGER   TO  4- 

INCH 

TOP. 

BD.  FT. 

-  TREES  8.0 

INCHES  D. 

B.H.  ANO 

LARGER  TO 

6-INCH  TOP. 
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YIELDS   PER   ACRE  OF  MANAGED,    EVEN-AGED  STANDS  OF   ENGELMANN   SPRUCE   AND   SUBALPINE  FIR 
SITE   INDEX     80,      30-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -      120.,    SUBSEOUENT  -  80. 
ENTIRE   STAND   BEFORE   AND  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 

AGE 
YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE  VOLUME 
CU.FT . 

SAWT I MBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

30 
30 

850 
505 

94 

72 

4.5 
5.1 

28 
29 

1010 
870 

340 
340 

0 
0 

345 

22 

140 

0 

0 

40 

498 

108 

6.3 

37 

1620 

1050 

0 

50 

485 

141 

7.3 

45 

2540 

1940 

0 

60 
60 

466 
172 

175 
78 

8.3 
9.1 

52 
52 

3660 
1740 

3040 
1510 

8700 
4800 

294 

97 

1920 

1530 

3900 

70 

167 

99 

10.4 

59 

2470 

2230 

8400 

80 

167 

123 

11.6 

64 

3350 

3100 

12800 

90 
90 

167 
61 

149 
80 

12.8 
15.5 

69 
70 

4390 
2420 

4120 
2320 

18200 
10800 

106 

69 

1970 

1800 

7400 

IOC 

61 

96 

17.0 

75 

3100 

2990 

14700 

110 

61 

113 

18.4 

79 

3660 

3550 

18100 

120 
120 

6 1 

20 

130 

63 

19.8 
24.0 

82 
83 

4320 
2060 

4210 
2040 

22200 
11200 

41 

67 

2260 

2170 

11000 

130 

20 

73 

25.9 

86 

2460 

2430 

13700 

140 

20 

84 

27.8 

89 

2890 

2860 

16600 

150 

20 

96 

29.7 

91 

3360 

3320 

19800 

TOTAL 

YIELDS 

9650 

8820 

42100 

MINIMUM 

CUTS 

FOR   INCLUSION  IN 

TOTAL   YIELDS — 

400.  CUBIC 

FEET 

AND  2000 

BOARD 

FEET 

MERCH . 

CU.  FT 

.  -  TREES 

5.0   INCHES  D.B.H. 

AND  LARGER   TO  4- 

INCH 

TOP  . 

BO.   FT.   -   TREES  8.0   INCHES   D.B.H.   AND  LARGER  TO   6-INCH  TOP. 

YIELDS   PER  ACRE  OF   MANAGED,    EVEN-AGED  STANDS   OF   ENGELMANN   SPRUCE   AND   SUBALPINE  FIR 
SITE    INDEX     80,      30-YEAR   CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL  -     120.,    SUBSEQUENT  -  120. 
ENTIRE   STAND   BEFORE   AND  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

30 
30 

850 
505 

94 
72 

4.5 
5.1 

28 
29 

1010 
870 

340 
340 

0 
0 

345 

22 

140 

0 

0 

40 

498 

108 

6.3 

37 

1620 

1050 

0 

50 

485 

141 

7.3 

45 

2540 

1940 

0 

60 
60 

466 
259 

175 
117 

8.3 
9.1 

52 
52 

3660 
2570 

3040 
2230 

8700 
7300 

207 

58 

1090 

810 

1400 

70 

249 

141 

10.2 

59 

3480 

3130 

11800 

80 

249 

173 

11.3 

64 

4680 

4300 

17800 

90 
90 

249 
103 

20S 
120 

12.3 
14.6 

69 
70 

5960 
3600 

5570 
3430 

24400 
15900 

146 

85 

2360 

2140 

8500 

100 

103 

142 

15.9 

75 

4540 

4360 

20900 

110 

103 

166 

17.2 

7B 

5410 

5220 

26300 

120 
120 

103 
35 

190 
96 

18.4 
22.4 

82 
83 

6310 
3160 

6120 
3120 

31900 
16900 

68 

94 

3150 

3000 

15000 

130 

35 

110 

24.0 

86 

3710 

3670 

20400 

140 

35 

125 

25.6 

88 

4300 

4260 

24300 

150 

35 

141 

27.2 

91 

4940 

4890 

TOTAL 

28700 
YIELDS 

11680 

10840 

52200 

MINIMUM 

CUTS 

FOR  INCLUSION  IN 

TOTAL  YIELDS- 

400.  CUBIC  FEET 

AND  2000. 

BOARD 

FEET 

MERCH. 

CU.  FT 

.  -  TREES 

5.0   INCHES  D.B.H 

ANO  LARGER   TO  4-INCH 

TOP. 

BO.  FT.  -  TREES  8.0  INCHES  O.B.H.   AND  LARGER  TO  6-INCH  TOP. 
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Three-Cut  Shelterwood 

YIELDS   PER  ACRE  OF   MANAGED,   EVEN-AGED  STANDS  OF   ENGELMANN   SPRUCE   AND  SUBALPINE  FIR 
SITE    INDEX     80,      30-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -        80.,    SUBSEOUENT   -  80. 
ENTIRE   STAND   BEFORE   ANO  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT, 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BO.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT  . 

TOTAL 
VOLUME 
CU.FT . 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BD.FT. 

30 

30 

850 
31* 

9* 

52 

*.5 
5.5 

28 
29 

1010 
680 

3*0 
3*0 

0 
0 

536 

*2 

330 

0 

0 

*0 

311 

78 

6.8 

38 

1260 

900 

0 

50 

305 

106 

8.0 

*5 

2030 

1650 

*200 

60 
60 

295 
1*9 

133 

80 

9.1 
9.9 

53 
53 

2920 
1810 

2530 
1610 

8300 
5700 

1*6 

53 

1110 

920 

2600 

70 

1*5 

99 

11.2 

59 

2510 

2310 

9200 

80 

1*5 

122 

12.* 

65 

3370 

3150 

13600 

<)0 
90 

1*5 
*0 

1*6 
6* 

13.6 
17.1 

70 
71 

*350 
1970 

*120 
1900 

18700 
9200 

105 

82 

2380 

2220 

9500 

100 

*0 

76 

18.7 

76 

2*10 

23*0 

11800 

110 

*0 

90 

20.3 

80 

2910 

28*0 

1*900 

120 
120 

*0 
13 

105 
*9 

21.9 
26.* 

83 
8* 

3*60 
1620 

3*00 
1610 

18*00 
9000 

27 

56 

18*0 

1790 

9*00 

130 

13 

57 

28.* 

87 

1930 

1910 

11000 

140 

13 

66 

30.* 

90 

2250 

2230 

13200 

150 

13 

75 

32.5 

92 

2620 

2600 

15700 

TOTAL 

YIELDS 

8280 

7530 

37200 

MINIMUM 

CUTS 

FOR   INCLUSION  IN 

TOTAL   YIELDS —     *00.  CUBIC 

FEET 

AND  2000. 

BOARD 

FEET 

MERCH. 

CU.  FT 

.   -  TREES  5.0  INCHES  D.B.H. 

AND  LARGER   TO  *- 

INCH 

TOP  . 

BD.  FT. 

-  TREES  8.0 

INCHES  D.f 

I.H.   AND  LARGER  TO 

6-INCH  TOP • 

YIELDS   PER  ACRE  OF  MANAGED,    EVEN-AGED   STANDS  OF   ENGELMANN   SPRUCE   ANO   SUBALPINE  FIR 
SITE   INDEX     80,     30-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -       80.,    SUB  SEQUENT  -  120. 
ENTIRE   STAND   BEFORE   AND  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

30 

30 

850 
31* 

9* 

52 

*.5 
5.5 

28 
29 

1010 
680 

3*0 
3*0 

0 
0 

536 

*2 

330 

0 

0 

40 

311 

7B 

6.8 

38 

1260 

900 

0 

50 

305 

106 

8.0 

*5 

2030 

1650 

*200 

60 
60 

295 
228 

133 
119 

9.1 
9.8 

53 
53 

2920 
2660 

2530 
2360 

8500 
8  500 

67 

14 

260 

170 

0 

70 

218 

1*1 

10.9 

59 

3520 

3220 

12800 

80 

218 

171 

12.0 

65 

*670 

*350 

18700 

90 
90 

218 
66 

20* 
96 

13.1 
16.3 

70 
71 

5990 
2930 

56*0 
2820 

25500 
13*00 

152 

108 

3060 

2820 

12100 

100 

66 

113 

17.7 

76 

3570 

3*50 

17300 

110 

66 

131 

19.1 

79 

*260 

*l*0 

21500 

120 
120 

66 
23 

151 
76 

20.  S 
2*. 6 

83 
8* 

5010 
2500 

*900 
2*80 

26300 
13700 

43 

75 

2510 

2*20 

12600 

130 

23 

B7 

26.* 

87 

2950 

2920 

16700 

1*0 

23 

100 

28.2 

89 

3**0 

3*00 

19900 

150 

23 

113 

30.0 

92 

3960 

3920 

23500 

TOTAL 

YIELDS 

10120 

9160 

*8200 

MINIMUM 

CUTS 

FOR   INCLUSION  IN 

TOTAL   YIELDS — 

*00.  CUBIC 

FEET 

AND  2000. 

.  BOARD 

FEET 

MERCH. 

CU.   FT.   -  TREES 

5.0   INCHES  D.B.H. 

AND  LARGER   TO  *- 

INCH 

TOP. 

80.    FT.   -   TREES   8.0   INCHES   D.B.H.    AND  LARGER  TO  6-INCH  TOP. 
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YIELDS   PER  ACRE   OF   MANAGED,    EVEN-AGED  STANOS   OF   ENGELHANN   SPRUCE   AND   SUBALPINE  FIR 
SITE    INDEX     80.      3C-YEAR   CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -     120..    SUBSEQUENT  -  80. 
ENTIRE   STAND   BEFORE   AND  AFTER   THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HE IGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE  VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BO. FT. 

30 
30 

850 
505 

94 
72 

4.5 
5.1 

28 
29 

1010 
870 

340 
340 

0 
0 

345 

22 

140 

0 

0 

40 

498 

108 

6.3 

37 

1620 

1050 

0 

50 

485 

141 

7.3 

45 

2540 

1940 

0 

60 
60 

466 
172 

175 

7B 

8.3 
9.  1 

52 
52 

3660 
1740 

3040 
1510 

8700 
4800 

294 

97 

1920 

1530 

3900 

70 

167 

99 

10.4 

59 

2470 

2230 

8400 

80 

167 

123 

11.6 

64 

3350 

3100 

12800 

90 
90 

167 

45 

149 
64 

12.8 
16.2 

69 
71 

4390 
1970 

4120 
1890 

18200 
8900 

122 

85 

2420 

2230 

9300 

100 

45 

78 

17.8 

75 

2460 

2  390 

11800 

110 

45 

92 

19.4 

79 

2990 

2910 

15100 

120 
120 

45 
15 

107 

53 

20.9 
25.4 

82 
84 

3550 
1730 

3480 
1720 

18500 
9500 

30 

54 

1820 

1760 

9000 

130 

1  5 

61 

27.4 

87 

2060 

2040 

11700 

140 

15 

71 

29.4 

89 

2430 

2400 

14100 

150 

15 

Bl 

31.4 

91 

2820 

2790 

16800 

TOTAL 

YIELDS 

9120 

8310 

39000 

MINIMUM 

CUTS 

FOR   INCLUSION  IN 

TOTAL  YIELDS- 

400.  CUBIC 

FEET 

AND  2000. 

BOARD 

FEET 

MERCH. 

CU.  FT 

.  -  TREES 

5.0  INCHES  D.B.H 

.   AND  LARGER   TO  4- 

INCH 

TOP. 

BD.   FT.   -   TREES   8.0   INCHES  D.B.H.   AND  LARGER  TO  6-INCH  TOP. 

YIELDS   PER  ACRE   OF   MANAGED.   EVEN-AGED  STANDS   OF   ENGELMANN   SPRUCE   AND  SUBALPINE  FIR 
SITE    INDEX     80,      30-YEAR   CUTTING  CYCLE 
THINNING  LEVELS=    INITIAL   -      120.,    SUBSEQUENT   -  120. 
ENTIRE   STAND  BEFORE   AND  AFTER   THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARSI 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWT IMBER 
VOLUME 
BO. FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

30 
30 

850 
505 

94 
72 

4.5 
5.1 

28 

29 

1010 
870 

340 
340 

0 
0 

345 

22 

140 

0 

0 

40 

498 

108 

6.3 

37 

1620 

1050 

0 

50 

485 

141 

7.3 

45 

2540 

1940 

0 

60 
60 

466 

259 

175 
117 

8.3 
9.1 

52 
52 

3660 
2570 

3040 
2230 

8700 
7  300 

207 

58 

1090 

810 

1400 

70 

249 

141 

10.2 

59 

3480 

3130 

11800 

80 

249 

173 

11.3 

64 

4680 

4300 

17800 

90 
90 

249 
75 

205 
96 

12.3 
15.3 

69 
71 

5960 
2910 

5570 
2780 

24400 
13000 

174 

109 

3050 

2790 

11400 

100 

75 

114 

16.7 

75 

3690 

3550 

17400 

110 

75 

134 

1B.1 

79 

4360 

4230 

21600 

120 
120 

75 
25 

156 
76 

19.5 
23.6 

82 
83 

5160 
2510 

5030 
2480 

26500 
13500 

50 

80 

2650 

2550 

13000 

130 

25 

88 

25.4 

86 

2970 

2940 

16500 

140 

25 

101 

27.2 

89 

3470 

3430 

19900 

150 

25 

115 

29.0 

91 

4010 

3970 

TOTAL 

23600 
YIELDS 

10940 

10120 

48000 

MINIMUM 

CUTS 

FOR   INCLUSION  IN 

TOTAL   YIELDS — 

400.   CUBIC  FEET 

AND  2000. 

BOARD 

FEET 

MERCH. 

CU.  FT 

.   -  TREES 

5.0   INCHES  D.B.H. 

AND  LARGER   TO   4- INCH 

TOP . 

BO.    FT.   -   TREES   8.0   INCHES  D.B.H.   AND  LARGER  TO  6-INCH   TOP . 
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Appendix  III:  Listing  of  Subroutine  ESSF 


SUBROUTINE  ESSF 

C 

C  LOCATION  FOR  ALL  SPECIES  -  SPECIFIC  STATEMENTS  APPLICABLE  TO  ENGELMANN 
C  SPRUCE  -   SUBALPINE   FIR    IN  COLORADO  AND  WYOMING. 

COMMON  ADD,AGE(2),AGE0,8A(2),BAS(2).8AS0,BAST,3AUS,BFMRCH,BFV0L, 
1CFV0L.0ATEI6I ,DBH(2)  ,  DBHE  ,  OBHO,  OBHT  ,  DEN  I  2  I  ,  DENO ,  DENT  ,  OMUS  ,  FB  A  (  2  )  , 
2FCTRI2 ) ,FDM(2 ) ,FDN{ 2) .FHT (2  I , FORETI 19) ,FVL( 2) ,HT( 2 ) , HTCUM , HTSO , 
3HTST,KAK,KN0,M[N, MNK, NBK.NC MP, NSUB.NWGP, PDBHE, PRET,PROD(2), REST, 
4SAVE,SBARB,SBARE,SBARG,SBAS,SITE,SLAND,TBA(2),TDM(2),TEM,TIME,  TM8R 
5,TMPO,TOT (2) ,TOTO.TOTT,TVL( 2 ) ,VOM( 2) . VLUS, OMR (  2) 

COMMON  ABFAG(5,15l,ACINTI5>,AOJ(5),AGETH(5,14),ALLCF(5,14l,ALOWC(5 
1  I  ,ALW8F I  5  I ,AMCAG(5,15  > ,A  MCUT(5, 14  1 ,AREA( 5,1*) , SOMA  I ( 5  I , BF AGE ( 5,15) 
2 ,BFINT 15  )  .CFAGEI5 , 15)  ,CFBF(  5,  14) .COMBFI 5  I  ,COMCU(  5  )  ,CUCY(  5  ),CUINT(  5 
3)  .CUM A  I ( 5  1 ,DBHTH(5, 14 ). DELAY! 5) , DENTH ( 5 , 14 1 , DL  E  V  I  5 ) , F  NBD I  5 ) , 
4FNCUI5I .GROWB (5 ,2 ,14) , GROWC ( 5 , 2 , 14 ) , GVL BF ( 5 ) , GVLCU ( 5),INVL(5,3,14) 
5  ,  NS I  15) , CPBD ( 5  > ,CPCU( 5),PAIBDI5),PAICU(5 ) ,POOR( 5 ) , R  EGN (5,3,14), 
6RINTI5),SARSP(5l,SeFI5),SHFLT(5,2,14), SHWOI 5, 2, 14  I , SMC ( 5 ) , SM S P ( 5 ) , 
7SUB8FI5,14),SUBCF(5,14) , SUMCF I  5 ) , S YST ( 5 ) , TH IN ( 5 ) , VLL V  I  5 , 3 , 14 ) , 
8WGNUM(  5  ),WGPDES(5,20),WGPN"M5,3)  ,  SPNUMI  5  ),TP8I5,7),PASP(5,7) 

COMMON  AC  BAR! 7) ,ARBK( 7) , B AR S I ( 7 , 14 ) , BFTH ( 7 , 27 ) ,CMTH(7,27) ,CUTA( 7,2 
17) ,CUTB (7 ,27)  ,HE LP ( 7 , 27 ) , NS BK ( 7 ) , OPEN ( 7, 27 ) , P BR S I ( 7 , 14 ) , PDCFN I  7 , 27 
2I,PDCFR(7,27),PSPLT(7,27),PUNC(7,27),SARETY(7,3  5),SLVG(7,27),SPLT( 
37,27),TMTY(7),UNCMLI7,27),PABR(7),PARTY(7,35I 

COMMON  ACFNL(5,7,15 ) , ACRGNI 5,7, 15 ) , ACS  1 1  5, 7, 14 ) , AC SP ( 5, 7 ) , GR BO ( 5 , 7 

COMMON  ACFNL(5,7,15) , ACRGN (5,7,15) , ACS I<  5 , 7 , 14 ) , AC SP ( 5 , 7 ) ,GRBO( 5,7 
1  ,15)  ,GRMC (5,7 ,15 ) , PS ( 5,7, 14  I ,STYP( 35 ) ,TYPNM( 35, 5 ) , PA S I ( 5 , 7, 14 ) 


COMMON   /BLKD/    IJ, IK,K I ,VOL,TVOL 

GO  TO    (10  ,20,30,40,50,60,70,80,90,  100,110,120  1,    I J 
SECTION   1  -   FIND  TOTAL  CUBIC   FOOT  VOLUME . 

10  D2H  =  DBH(IK)    *   DBH(IK)    *   HT ( IK) 
IFID2H   .GE.    22500. C)    GO  TO  11 

TOT(IK)  =  (0.0023575  *  D2H  -  0.00524  *  BAS(IK)  +  0.34480)  *DEN(IKI 
GO  TO  12 

11  TGT(IK)    =    (0.0020103   *  D2H  -  0.01052   »  BAS(IK)    ♦   7.60196)*  DEN(IK) 

12  RETURN 

SECTION  2  -  VOLUME  CONVERSION  FACTORS. 

KERCH.   CU.    FT.   -  TREES  6.0   INCHES   D.B.H.   AND  LARGER   TO  4- INCH  TOP. 
SO.   FT.   -  TREES   10.0   INCHES   D.B.H.    AND  LARGER   TO   8-INCH  TOP. 

2C  00  21   J  =  l  ,2 

FCTR(J)  =  0.0 
21  PROO(J)    =  0.0 

00  26   1=1, KNO 

IF(VDMII)    .LE.   4.991    GO  TO  26 

FCTRIII    =   0.82375   +   3.45569   /   VDM(I)    +   0.C0313   «   VDM(I)    *  VDM(I) 
1  -   28.86783   /    IVDM(I)    *  VDM(Il) 
IFIFCTR(I)    .GT.   G.99)    FCTRIII   =  C.99 

24  IF(VDMII)    .LE.   7.99)    GO  TO  26 
IF(VDM( I )    .GE •    16.5)    GO  TO  25 

PROD(I)    =  4.59159  -   214.06370   /    (VDM(I)    *  VDMID)    *■  0.39782  * 

I  ALOG10(8A( I ) ) 
GO  TO  26 

25  PROOI I  I    =  8.59422 

II  )  ) 

26  CONTINUE 
RETURN 


19.545C7   /   SQRTIVOM(I))    t   0.44432   *  AL0G10IBA1 


C  SECTION  3  -  GROWTH  FOR  NEXT  PERIOD. 
C 

30  DO  35  1=1,2 

TMDY  =   AGE ( I )    +  TIME 
IFITMDY   .LT.    TEM)   GO   TO  35 

FOMII)  =  1.62917  +  1.03371  *  DBH(I)  *  0.01304  *  SITE  -  0.90669  * 
1  ALOGIO(SBAS) 

31  IF(DBHII)    .GE.   10.0)   GO  TO  32 

FDN(I)    =    (-0.03967   +   0.CO382   *   DBH1I)    *   SBAS)   *  0.01 

IF(FDNII)    .LT.   0.0)    FDN(I)    =  O.C 

FON ( I )    =   DEN ( I )    *    (1.0  -   FDN ( I ) ) 

MNK  =  FDN  I  I )    +  0.5 

FDN ( I )    =  MNK 

GO  TO  33 

32  FDN I  I )    =  DEN ( I ) 

33  F8AII)    =   0.0054542   *   FDM(I)    *   FDM(I)   *   FDN ( I ) 
IF(AGEII)    .GE.    100.01    GO  TO  37 

FHT ( I )    =  -13.71751   +   0.15087   *   SITE   +   0.00126  *   AGE ( I )   *   AGE ( I )  * 
1  0.01371   *   AGE ( I )    *   SITE   -  0.000C6   *   AGE  I  I )    *   AGE (  I  )   *  SITE 
GO  TO  38 

37  FHT ( I )    =   0.91859  -   100.43601   /   AGE ( I )    +  0.62318   *   ALOGIO(SITE)  + 
1   40.08154  *   ALOGIO(SITE)    /   AGE ( I ) 

FHT ( I )    =   10.0  **  FHT ( I ) 

38  D2H  =  FDM(I)    *  FDM(I)    *  FHT ( I ) 
IFI02H   .GE.   22500.0)    GO  TO  34 

FVLII)  =  (0.CO23575  *  D2H  -  0.00524  *  FBA(I)  +  0.34480)  *  FON(I) 
GO  TO  35 

34  FVL(I)    =    (0.0020103   «D2H  -  0.01052   «   FBA(I)    +   7.601961   *  FDN(I) 

35  CONTINUE 
RETURN 

C 

t  SECTION  4  -   FUTl'RE  UNTHINNED  UNOERSTORY   IF   OVERSTORY   REDUCED  NOW. 

40  DMUS  =   1.62917  +   1.03371   *  0BHI2)    +  0.01304  *   SITE  -   0.90669  * 
1  AL0G10(BAS(2)> 


IFIDBHI2)    .GE.   10.01    GO  TD  41 

ONUS   =    (-0.03967  +   0.00382   »   0BHI2)    *   BAS ( 2 ) )    *  0.01 

IFIDNUS   .LT.   0.0)   ONUS  =  0.0 

ONUS  =   DEM2)    *    (1.0   -   DNUS  ) 

MNK   =   DNUS  ♦  0.5 

DNUS  =  MNK 

GO  TO  42 

41  ONUS   =  0ENI2) 

42  BAUS  =   0.0054542   *  OMUS   *   DMUS   *  DNUS 
IFIAGEI2I    .GE.    100.0)   GO  TO  47 

HTUS  =  -13.71751   ♦  0.15087   *   SITE   ♦  0.00126  *   4GEI2)   *   AGE ( 2 )  ♦ 
1   0.01371   *  AGE ( 2 )    *   SITE   -   0.030C6  «■   AGE ( 2 )    *   AGE ( 2 1    »  SITE 
GO  TO  48 

47  HTUS   =  0.91859  -   ICO. 43601   /   AGE ( 2 )    ♦  0.62318   *   AL0G13ISITE)  ♦ 
I   40.C8154  *   AL0G10ISITE)    /   AGE ( 2 ) 

HTUS   =   10.0  »*  HTUS 

48  D2H   =   DMUS   *   DMUS   *  HTUS 
IFID2H  .GE.   22500.0)   GO  TO  43 

VLUS  =    (0.0023575   *   02H  -  0.00524   *   BAUS   +   0.34480)   *  DNUS 
GO  TO  44 

43  VLUS   =    (0.C020103  *  D2H  -  0.01052   »   BAUS  +   7.60196]    *  ONUS 

44  RETURN 

SECTION   5   -   NEW  D.B.H.    AFTER  THINNING. 

50  IFIPRET   .LT.    50.0)    GO  TO  51 

DBHE    =   0.02666  +   1.33655   *   OBHO  - 

GO  TC  52 

51  PDBHE   =   0.33206  *   0.98346   *   ALOGIO(DBHO)   -  C. 14710 
DBHE   =   10.0  **  PDBHE 

52  RETURN 


00306   *   PR  ET    *  DBHG 

ALOGIO(PRET) 


C   SECTION   6  -  CUBIC   FEET   AS   BYPRODUCT   OF   SAWLOG  CUT. 

C 

60  ADD   =   VCL   *    (1.54375   +   11  .43324   /  PPHIKI)) 
ADD  =   TVCL  -  ACD 

I F ( ADD   .LT.   COMCU(KAKI)    A  DO   =  0.0 
RETURN 

C    SECTION   7  -   VOLUME    IF   THINNED  NOW   AND    IF   THINNED   IN   TIME  YEARS. 

7:   HT(KI)    =   HT(KI)    -   1  .8  1392   +   1  7.77922   *   SQRTd.O   /   PR  E  T  ) 
TEM   =   TBA(IK)    /    (0.0054542   *   TDM ( IK )    *  TDM(lKl) 
D2H   =   TDM(IK)    *   TDMIIKI*  HT(KI) 
IFID2H   .GE.    22500.")    GO  TO  71 

TVL(IK)    =    (0.0023575   *   02H  -   0.00524   »  TBA(IK)    +   0.34480)   *  TEM 

GO  TC  72 

71  TVL(IK)    =    (C. 0020103   *  D2H  -   0.01C52   *  T8AIIK)    +   7.601961    *  TEM 

72  RETURN 


C   SECTION   8  -   STATUS  AT   END  OF   PERIOD   IF   THINNED   AT   START  OF  PERIOD. 

C 

SO   J   =   TIME   /  RINT(KAK) 
DO  33    1=1, J 

IFIT3AI1)    .LE.   0.0)    GC  TO  83 

HTll)    =   HTIll   -   1.81392   +   17.77922   *   SORT (1.0   /  SAVE) 
FDM(i)    =   1.62917   +   1.03371   *   TDM(l)    +   0.01304   *   SITE   -   0.90669  * 
1   ALOGICITBA(I) ) 
I F ( AGE ( K 1 )    .GE.    10C.0)   GO   TO  87 

FHT(l)    =  -13.71751    ♦   0.15087   *   SITE   *   0.00126   *   AGE(KI)   *   AGE ( K I ) 
1    *   0. 01371   *   AGE ( K I )    *  SITE   -  0.00006   *   AGE ( K I )    *   AGE ( K I )   *  SITE 

GO  TC  38 

87   FHT ( 1 )    =   0.91859  -   100.43601    /   AGE(KI)    +  0.62318   *  AL0G10ISITE)  ♦ 
1    40.08154  *  ALOGIO(SITE)    /   AGE ( K I ) 

FHT ( 1 )    =   10.0   **   FHT ( 1 ) 
BS   MNK   =    (TBA(l)    /    (0.0054542   *   TDM ( 1 )   *   TDM ( 1 ) ) )    +  0.5 

IF(TDMIl)    .LT.    10. C)    GO   TO  81 

FDN 1 1 )    =  MNK 

GO   TO  82 

91   FDN(l)    =    (-0.03967  +   0.00382   *  TDM(l)    *   TBA(D)   *  0.01 
IF(FDNd)    .LT.   0.0)    FCN(l)    =  0.0 

TFM  =  MNK 

NKM   =   TEM   *    (1.0   -   FCM11  I    *  0.5 
FDN ( 1 )    =  NKM 

82   FBAI1I    =   FON(l)    *   0.0054542   *   FDM(l)    *   FDM( 1 ) 
TDM ( 1 )    =  FDM(l) 
TBI ( 1 )    =  FBA(l) 
HTll)    =   FHT ( 1 ) 

AGE ( K I )    =   AGE  I K I )    +  RINTIKAK) 
S3  CONTINUE 

D2H  =  FDMI1)  *  FDMI1)  *  FHT ( 1 ) 
IFID2H   .GE.   22500.0)    GO  TO  84 

FVLII)    =    10.0023575   *   D2H  -  0.00524  *  FBA(l) 

GO   TO  85 

84  FVL(l)    =    (0.0020103   *   D2H   -  0.01C52   *  FBA(l) 

85  RETURN 


0.34480)  *  FDN ( 1  ) 
7.60196)   *   FDN ( 1 ) 


SECTION   9  -  HEIGHT   AND  VOLUME   BEFORE  THINNING. 

90  I F ( AGE C   .GE.    100.01    GC  TO  91 

HTSO   =  -13.71751   +  0.15087   4   SITE   +   0.00126   *   AGEO  *   AGEO  + 
1   0.01371   *   AGEO  *   SITE   -  0.00006   *   AGEO  *   AGEO   *  SITE 
GO  TO  92 

91  HTSO  =   0.91859  -   100.43601   /   AGEO   +  0.62318   *   ALOGIO(SITE)  ♦ 
1   4Q.G8154   *   ALOGIO(SITE)    /  AGEO 

HTSO  =  10.0  »*  HTSO 
9?   HTSO  =   HTSC  *  HTCUM 

D2H   =  OBHO  *   DBHO  *  HTSO 
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IFID2H   . GE .   22500.0)    GO  TO  93 

TOTO  =    (0.0023575   *  02H  -  0.00524   *   BASO   +  0.344801   *  OENO 
GO  TO  94 

93  TOTO  =    10.0020103   *  02H  -   0.01052   *  BASO  ♦   7.601961    *  DEMO 

94  RETURN 

SECTION   10   -  HEIGHT   AND  TOTAL   CUBIC   FEET  PER   ACRE   AFTER  THINNING. 

100  AOOHT   =  -   1.81392  +   17.77922   *   SQRTd.O  /   PRET ) 
HTCUM   =  HTCUM  +  ADDHT 

HTST   =   HTSO  ♦  ADDHT 

D2H   =  DBHT  *  DBHT   *  HTST 

IFI02H   .GE.  22500.01   GO  TO  101 

TOTT   =    (0.0023575   *  D2H  -  0.00524  *  BAST   ♦   0.34480)   *  DENT 
GO  TO  102 

101  TOTT   =    10.0020103   «  D2H  -  0.01052   *   BAST   t   7.60196)   *  DENT 

102  RETURN 


SECTION   11   -   D.8.H.    AT   END  OF   PROJECTION  PERIOD. 

110  DBHO  =   1.62917  +   1.03371   *   DBHT   +  0.01304  *   SITE  -   0.90669  * 
1   AL0G10 ( BAST  > 
RETURN 

SECTION   12   -   MORTALITY   AS   A  PERCENTAGE  OF   INITIAL  DENSITY. 


120   DENO  =    (-0.03967   +  0.00382   *   DBHT   *   BAST)   *  0.01 
RETURN 
END 
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